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Abstract 

We consider the 2x2 parabolic systems 

on a domain {t, x) £]0, +oo[x]0, /[ with Dirichlet boundary conditions imposed at a; = and 
at x = The matrix A is assumed to be in triangular form and strictly hyperbolic, and the 
boundary is not characteristic, i.e. the eigenvalues of A are different from 0. 

We show that, if the initial and boundary data have sufficiently small total variation, then 
the solution u'^ exists for all f > and depends Lipschitz continuously in on the initial and 
boundary data. 

Moreover, as £ — > 0"*", the solutions u'^it) converge in to a unique limit u{t), which can 
be seen as the vanishing viscosity solution of the quasilinear hyperbolic system 

ut+ A{u)uj,=Q, xe]0, 

This solution u{t) depends Lipschitz continuously in w.r.t the initial and boundary data. We 
also characterize precisely in which sense the boundary data are assumed by the solution of the 
hyperbolic system. 
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1 Introduction 

This paper deals with the initial-two-boundaries value problem 

Ut + A{u)ux = 0, xe]0, te]0, -|-oo[ 

M(0,a;) = uo(a;), (1.1) 

M(t, 0) = Uho(^), u{t,l) ^ uu{t). 

The crucial hypotheses we assume are that the matrix A is strictly hyperbolic with eigenvalues 
different from and that the initial and boundary data are small in BV norm and close to a 
constant state u* . 

An existence result for hyperbolic boundary value problems was proved in [JSJ |2] using an 
adaptation of the Glimm scheme introduced in |23 ■ Improvements of the results in j2Sl E] have 
been obtained by a wave- front tracking technique introduced in [2j and later used in a series of papers 
(^lElElEIEl^SEl) to establish the well posedness of the Cauchy problem. Such a wave-front 
tracking technique was adapted to the initial-boundary value problem in PP, where a substantial 
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improvement of the results in |25[l31j was achieved. The weU posedness of the initial-boundary value 
problem was then proved in j20| relying on the wave- front tracking technique described in pj. 

All the results quoted so far deal with conservative systems; a comprehensive account of the 
stability and uniqueness results for the Cauchy problem for a system of conservation laws can be 
found in We refer, instead, to and to for a general introduction to the systems of 

conservation laws. 

In 0] 121 El and (71 a different problem was dealt with: let be a family of solutions to the 
parabolic systems 

ul + A(u'')ul = £<^. 

One expects that as e ^ 0+ the solution u'^ converges in some sense to a solution of the corresponding 
hyperbolic system 

ut + A{u)ux — 0. 

The mathematical proof of this convergence was obtained via a suitable decomposition of the gradient 
of the solution along travelling waves. We refer to (2j for an account of the proof of the convergence 
of the vanishing viscosity approximation and of the uniqueness and the stability of the vanishing 
viscosity limit: it is important to underline, however, that in 7* the systems considered are not 
necessarily conservative. 

The vanishing viscosity approximation of initial-boundary value problems was studied in numer- 
ous works: in the following, we will briefly refer to some of the principal results, without any sake 
of completeness. Moreover, if not otherwise stated, the systems considered are supposed to be in 
conservation form. 

In particular, in '32' it was considered the vanishing viscosity approximation 

of an initial-boundary value problem and it was given a precise description of the first term of the 
expansion of in the neighborhood of a point where two shocks or a shock and a boundary layer 
profile meet. 

The works |2H (23 dealt with the general parabolic approximation 

ul + f{u')x^e{B{u')ul)^, (1.2) 

where the viscosity B{u) is invertible but in general different from the identity. It was proved the 
existence of a T > such that u' converges in L°°((0, T); L^(R+)) to a solution of 

Ut + .f{u)x = 

and it is given a precise characterization of the boundary condition induced in the hyperbolic limit. 

In p4) it was introduced an Evans function machinery to study the stability of boundary layer 
profiles: the parabolic approximation considered was in the form H1.2|l , in the case of an invertible 
viscosity matrix B and of a non characteristic boundary (i.e. all the eigenvalues of Df{u) were 
supposed to be different from zero). However, the analysis was extended in a series of paper (|3S1 
EHlEniEni) to the boundary characteristic case and to very general parabolic approximations, with 
non invertible viscosity matrices. 

In (31 it was considered the family of initial-one-boundary value problems 

{uf + A{u^)Ux = eu%x, a; e ]0, -l-oo[, te]0, -|-oo[ 
U''{0,X) = Uq{x), U^{t,0) = Ub{t), 

it is proved the (global in time) convergence of approximated solutions and the stability and the 
uniqueness of the limit. In 0| the boundary characteristic case was allowed (i.e. one characteristic 
field was allowed to have speed close to that of the boundary) and the crucial tool in the proof of 
the convergence and the stability is the introduction of a suitable decomposition of the gradient of 
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the vanishing viscosity solution. Moreover, we underline that, as in 0, the systems considered were 
not necessarily in conservation form. 

In the present paper we will consider the vanishing viscosity approximation for the initial-two- 
boundaries value problem: 

uf + A(u^)u% ^ eu%^, xe]0,l[, te]0,+oo[ 

u^O,x)^uo{x), (1.3) 

^ u'{t,0) = uboit), u'{t,l)^uu{t). 
We will assume that A is in triangular form, i.e. 

Ai(ui) 

A{u) =1 I , (1.4) 

g(ui, U2) A2(mi, U2) 

and sufficiently smooth in a compact neighborhood K oi a. fixed point u* . Moreover, we assume A 
to be uniformly strictly hyperbolic, in particular we assume that there exists a constant c > (2c is 
then the "separation speed") such that 

Ai(m) < -c < < c < A2(u) MueK. (1.5) 

The above condition means that the speed of the boundary (in our case 0) is strictly different from 
the characteristic speeds of the two families of waves. 

We denote with ri{u) the first eigenvector of A{u), corresponding to the eigenvalues Ai(u), and 
with r2 the second one. Due to the particular structure of A, we normalize ri and r2 as 

((l,0),ri(^.)) = l, ^2-(j). (1.6) 

The dual base of {ri{u), r2) is denoted by (£1, £2{u)). 

We will assume that the initial data uq and boundary data Ubo, Ubi have sufficiently small total 
variation, i.e. 

TotVar(?2o), TotVar(ufco), TotVar(ufc;) < Si (1.7) 

for a suitable Si << 1. Moreover, since we will study boundary layers with small total variation, we 
assume that there exists a value u* such that 

ll-So - U*||oo < (51 ||U60 - '"*||oo < (^i - U*|joo < ^1- (1-8) 

For technical reasons, we will also assume some stronger regularity: the boundary and initial data 
will be sufficiently smooth and will satisfy 

\\d^uo/dx^\Li, \\d^Ubo/dt^\\Li, Wd^Ubi/dt^Li < M < +x, j = 2, . . . n, (1.9) 

for some n G N and some large constant M . Some observations about the extension of our results 
to the case of boundary and initial data with weaker regularity will be made in R emark 1 1.21 

We will denote by Uq, Ub the set of functions uq, Ub satisfying (|1.7|l . H1.8|) . H1.9|l in ]0, Z[ or 
]0, +oo[, respectively. We also define the sets I?o Q i^(0, l)^ 'Db Q Ll^^{0, +oo) of functions such 
that 

TotVarjuo} < Si, TotVar{ufc} < Si, (1.10) 

respectively. 

Remark 1.1. The fact that we will consider only 2x2 triangular systems does not affect very deeply 
the structure of the problem, but leads to some considerable simplification in the computations. In 
particular, since the matrix A is in triangular form, we will see in Section |21 that the generalized 
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eigenvector of the travelling wave profile of the second family is constant, and so it is the generalized 
eigenvector of the boundary layer profile of the second family: such a feature simplifies the compu- 
tation of source terms, which is performed in the Appendix IA.2.11 Since also the expression itself of 
the source terms is simpler, the consequent estimates, carried on in Sectional are easier in the case 
of a triangular system than in the general one. 

We refer, instead, to Remark 1 1.21 for some considerations about the hypotheses of regularity we 
have assumed. 

The first theorem concerns the existence of a solution to the parabolic problem (|1.3|) ; moreover, 
it ensures that such a solution satisfies stability estimates independent on e. 

Theorem 1.1. Suppose uq € Uq, Ubo-, Ubi € Ub and A is of the form (|1.4() and satisfies (|1.5(l . 
Then, for any e > 0, the system (|1.3|l has a unique solution u^{t) defined for all t > 0. 

This solution depends Lipschitz continuously in on the initial and boundary data: indeed, let 
vq £ Uq, Vbo, Vbi £hlb be the initial and boundary data of a solution v^{t) of p. 3(1 . Then for some 
constants Li and L2, depending only on the matrix A and the bound on the initial and boundary 
data Si , the following holds: 



\\v%t) - U%t)\\Li < Li(^\\vo ~ UoIIlHO,;) + ll^fcO - Who|Ui(0,+oo) + Ikw - U6i||Li(0,+oo) 

+ L2f|<-s| + |Vi-%/i| 



(1.11) 



The second theorem concerns the limit as £ — > 0^. Since we have a uniform bound on the total 
variation, by Helly's theorem there is a subsequence of converging in to a limit function u(t') 
on a countable dense set of times t„. By the stability estimate Hl.ll|l . the convergence is on the 
whole M+. 

However, different subsequences could a priori converge to different limits: we will actually prove 
that the limit is unique and that moreover the semigroup property holds. 

Theorem 1.2. As e — > 0^, the sequence u'^it) of solutions of (|1.3|l converges to a unique function 
u{t) for all t > 0: we denote such a limit by 

u[t) =_Pf[?2o, Ubo, Ubi]- 

This convergence defines a unique semigroup 

S : [0, +00] X Uq X Kb ^ Ub Vq X Kb >^ Ub 

{t, Uo, Ubo, Ubi) [pt[uQ, Ubo, Ubi], Ubo{ - +t), Ubi{- 



which satisfies the following stability estimates in L^{0,1): 



Pt[Uo, Ubo, Ubl\ -ps[Vi3, Vbo, Vbl 



<Li \\vo~ -"01^1(0, 1) + \\vbO - ■"&o||l1(0, +00) 



+ \\Vbl - ww|Ili(o,+oo) I + L2\t 



(1.12) 



(1.13) 



for some constant Li, L2 depending only on A and on 5i. 

Remark 1.2. By the stability estimate H1.13|l the semigroup S defined by H1.12|l can be extended 
to initial and boundary data that satisfy much weaker regularity assumptions, i.e. uq e Vq and 
Ubo, Ubi G T>b- Indeed, let {pk\ be a sequence of regularizing kernels and let uq £ Vq. Then pk * mq, 
Pk*Ubo and pk*Ubi are initial and boundary data that satisfy the hypothesis l|1.9|l : they are smooth 
and 

\\d{pk * uo)/dx\\Li < Tot Var{t2o} < Si \\d{pk * Ubo)/dx\\Li < Si \\d{pk * Ubi)/dx\\Li < 5i 
d^pk * UQ)/dx^\\^, = d(^{d^~^pk/dx^~^) * UQ^/dx < M{k, j)Si j = l,...n 

d^ (^Pk * Ubo^/dx^ ^^<M{k,j)Si d^ (^pk * UQ)j / dx^ ^^<M{k,j)Si j = l,...n. 
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The last estimates ensures that, for any fixed k, the norm of the derivatives is finite: the bound 
is not uniform with respect to k but, since the constant Li in H1.13|l does not depend on the bound 
M in (|1.9I) . it is enough to prove the extendibihty of the semigroup to the whole domain Dq. Indeed, 
let u% the sequence of solutions to the systems 

< m|(0, x) = pk* uo 

ul{t, 0) = pk*Ubo ul{t, I) = pk*Ubi 

Theorem II .21 ensures that, for any k £ N and for any t > 0, the sequence converges as e ^ 0+ 

to some limit function we will call Uk{t). Then Uk{t) is a Cauchy sequence since by (|1.13ll 

\\uk{t)-Uhit)\\LHQj) < il(||(pfc-p/i)*Uo||Li(OJ) + ll(/'fc-P/i)*Wfto||Li(0,+oo) + ll(Pfc-p/0*^bniLMO,+oo))- 

The same estimate H1.13|l ensures that the limit limfc^+oo Uk(t) does not depend on the choice of 
the sequence pk and therefore the extension 

Pt[uo, Ubo, Ubi] = lim u''{t) 

k — * + CxD 

is well defined. 

For simplicity, in the following we won't prove that, if (uq, Ubo, Ubi) belongs to VoxiybX Vb but 
not to i^o X X Z/^b, then the solution of the system converges as £„ 0+ to pt[uo, WbOi Ubi]- 
However, we will exploit the extendibihty property described before, in particular in Section 16.11 
we will consider the vanishing viscosity solution of the Riemann and of the boundary Riemann 
problem, actually meaning the extension of the semigroup of the vanishing viscosity solution to 
piecewise constant initial and boundary data. 

The function u{t) — pt[uf), Ubo^ Ubi] is the vanishing viscosity solution to 

Ut+A{u)ux = 0. (1.14) 

Note that it is not a weak solution, unless the system is conservative, but one can prove that it is a 
viscosity solution, in the sense of ^2.. In particular, we obtain that, for a.e. t, the limits 

lim u{t, x) = u{t, 0+), lim u{t, x) = u{t, P) (1.15) 

and the boundary data Ubo(t), Ubi{t) can be connected by boundary profiles, i.e. there exists a 
solution of the boundary value problem 

f A{v)vx^Vxx, a;e]0, +oo[ f A{v)vx ^ v^^, x e] - oo, 0[ 

\ v{0) = Ubo{t), lim^j^+oo v{x) = u{t, 0+) \ v{0) = Ubi{t), lim2,^_oo v{x) = u{t, l^) 

respectively. This means that the boundary datum Ubo lies on the stable manifold of u{t, 0+), and 
the boundary datum Ubi lies on the unstable manifold of u{t, l^). 
The paper is organized as follows. 

First of all we make a change of variables in H1.3|l : let u{x, t) := u^{x/e, t/e). Then is 
equivalent to the system 

ut+ A{u)ux ^Uxx, a; G ]0, L[, t e ]0, +oo[ 
< u{{),x) =uq{x), (1.16) 
u{t,Q) = Ubi:){t), u{t, L) ^ UbL{t) 
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where L — l/e, Ubo{t) = UbQ{t/e), UbL{t) — Ubi{t/e), uo{x) — UQ{x/e). One can easily check that 

TotVarjubo} = TotVarjMbo} < Si TotVarjuM} = TotVar(wbL) < Si 
Tot Varjuo} = Tot Varjuo} < Si. 

Moreover, the derivatives of the boundary and initial data satisfy 

Wd^uo/dxmLi, Wd^Ubo/dt^h^, Wd^UbL/dtmLi <Ms^~' <Si j^2,...n (1.17) 

for £ small enough. 

The crucial tool in the proof of the convergence of the solution of as the scaling parameter 

e 0+ is Helly's theorem. One needs therefore to prove a uniform bound on the total variation, 
independent on the length of the interval L and on the norm of the boundary and initial data. 

In Section|21we prove a priori bounds on the solution of (|1.16|l that ensure the local existence and 
smoothness of solution. Moreover, we will show that, as long as the total variation of the solution 
remains small, the norm of u^x is small too and the solution itself can be prolonged in time. 
The proof is based on the following observation: (|l.lt)|) can be seen as a perturbed heat equation 
and therefore one is led to introduce suitable convolution kernels. Since the technique used in this 
section does not depend on the dimension of the solution u, we perform the computations for the 
n X n system. 

In SectionOlwe introduce the crucial tool in the proof of the BV estimates: a suitable decompo- 
sition of the gradient of the solution. In the boundary free case , the gradient Ux is decomposed 
along a suitable set of unit vectors fi, i = 1, . . . , n, which correspond to the tangent vectors of the 
travelling wave profiles of 

Ut + A{u)Ux = Uxx- 

In the single boundary case [2, instead, the gradient Ux is decomposed along n travelling wave 
profiles (the same as in the boundary free case) and along a boundary profile, i.e. a solution to the 
stationary system 

Uxx = A{u)ux. 

Such a boundary profile lays on a manifold whose dimension is related to the number of negative 
eigenvalues of A{u), i.e. to the number of characteristic fields that leave the domain x > 0. 

In our case, the basic idea is to split the part of the gradient due to the presence of the initial 
datum from the part due to the boundary data: the first part will be decomposed along the same 
tangent vectors fi, f2 to travelling wave profiles introduced in 7\. Moreover, following the same 
ideas as in in order to decompose the part of the gradient due to the boundary data we use 
double boundary profiles, i.e. suitable solutions of the stationary system 

i"^^4V^ (1.18) 
[ Px^ A{u)p. 

In the linear case the two components of the system H1.18|l are decoupled and one can show that 
there is a solution of the boundary value problem 

{Ux =P, 
Px = A{u)p, (1.19) 
uiO) = C/fco, u{L) = UbL 

with total variation uniformly bounded with respect to L. 

In the general case, the idea is to emulate the linear case, using the center-stable manifold theorem 
coupled with a contraction mapping argument: one finds that, provided the difference \UbQ — Ubhl 
is small, there is a solution of H1.19|l with uniformly bounded total variation. Such a solution can be 
seen as the sum of two components, one exponentially decreasing as a; — > +oo, the other as a; — > — oo: 
we will denote by f i and f2 the tangent vectors to the first and the second part respectively. It is 
important to underline, however, that in the non linear case the two components are coupled: indeed, 
one finds that Ai, the speed of exponential decay of the first component, depends also on the second 
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component, and viceversa A2 depends on the first component. The introduction of the generaUzed 
eigenvalues Ai and A2 allows the equations satisfied by the components of the decomposition to be 
exactly in conservation form. 

The decomposition of the gradient along travelling waves profiles and double boundary layers 
takes the form 

Ux ^ Vlh + V2f2 + Plh + P2f2- (1-20) 

In Section I3. II we will show that, because of the triangular structure of the matrix A, the vector f2 
and f2 can be chosen to be identically equal to r2 = (0, 1) and Ai is identically equal to Ai. 

Note that H1.2()|l is a system of 2 equations in 4 unknowns: this allows some freedom in choosing 
in the most suitable way the boundary and initial conditions. The precise expression of all the 
boundary and Cauchy data we will impose on vi, V2, pi and p2 can be found in Section [3. 31 in the 
following however we will sketch the crucial ideas involved in the choice of those conditions. 

We need a preliminary observation: besides that in the choice of the boundary conditions, some 
freedom is also allowed in the attribution of the source terms. Indeed, if one inserts (11.20(1 in the 
system 

Ut + A{u)Ux - Uxx = 

obtains the equations 

Vu + (Ait;i):r - Vixx + Pit + {>^lPl)x - Plxx = 

V2t + {>^2V2)x - V2xx + P2t + {>^2P2)x - P2xx = Sl{t, x) 

for some function si whose exact expression can be found in the ADDendix lA.2.1l and is not important 
at the moment: however, it is crucial to observe that it is identically zero when the solution is exactly 
a travelling wave or a double boundary profile. Moreover, in general such a source term is spread 
on the whole interval ]0, L[: since p2, the part of the double boundary layer exponentially decaying 
as a; ^ — oo, should be affected only by the datum in x — L, it seems reasonable to impose 

Vu + {XlVi)x - Vixx = Pit + {XlPl)x - Plxx = 

(1.21) 

V2t + iX2V2)x - V2xx = Slit, x) P2t + {X2P2)x - P2xx = 0. 

As regards the boundary and initial data we impose on the components pi, p2, vi and U2, we first 
observe that, since pi and p2 are the components of Ux along double boundary profiles, we don't 
want them to be influenced by the initial datum. Hence we impose 

pi(0, x) = P2(0, x) = 0. 

Moreover, pi is the exponential decreasing component of the boundary profile and hence it should 
not be affected too much by the datum on the boundary x = L: more precisely, since the goal is 
to establish a uniform bound on the norm of pi, it seems reasonable to look for some boundary 
condition that minimizes the increment of ||pi||li(o, l) due to the datum on the boundary x = L. 
An integration by parts ensures that 

— j \px{t, x)\ < \pix - XiPi\{t, L) + \pxx - Aipi|(t, 0) 
and therefore we will impose 

\pix - \iPi\{t, L) = 

and, by analogous considerations, 

\P2X - X2P2\it, 0)eeO. 

On the other hand, vi and V2 are the components of Ux along travelling profiles and therefore we 
don't want them to be strongly influenced by the presence of the boundary data. We observe that, 
in the hyperbolic limit 

Ut + A{u)ux = 0, 
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the waves of the first family go out from the domain through the boundary a; = 0: we would like to 
emulate such a behavior in the parabolic approximation. More precisely, since the aim is to show 
a uniform bound on the norm of vi, we look for some boundary condition that ensures that the 
derivative of the wave in the parabolic approximation crosses the boundary, as in the hyperbolic 
limit. To make the situation clearer, it is useful to consider the simple examples that follow: consider 
the linear scalar equation 

zt + Xlz^ - z^^ = (1.22) 
with some Dirichlet condition imposed on the boundaries a; = and x = L, for example 

z{t, 0) = 0, z{t, L) = 1. (1.23) 

Moreover, let z^{t, x) be a solution of H1.22|l and (|1.23|l : the initial condition is not important at 
the moment, but suppose for simplicity that Tot Var{z''^(0, x)} — 1. For sure Tot Var{2;^(t)} > 1 
and hence the derivative of z^ cannot cross the boundary a; = 0, or at least the loss of total variation 
that occurs at a; = has to be compensated by an increase at a: = L. 

On the other hand, let z^ {t, x) be a solution of p. 22(1 that satisfies a homogeneous Neumann 
condition at a; = 0, for example 

Zf (t, 0) = 0, Z^(t, L) EE 1, 

then an integration by parts ensures that 

|^%f(t, x)\dx<-\z^,At, 0)1, 

and hence the total variation of z^ is flowing out from the domain through the boundary a; = 0. 

Hence we are are led by the previous considerations to impose on the boundary x = a homo- 
geneous Dirichlet condition on the function wi, which corresponds to the derivative of a travelling 
wave of the first family: 

vi(t, 0) = 0. 
The considerations that motivate the choice 

V2{t, L) =0 

are completely analogous. 

In Section 01 we exploit the decomposition H1.2()|l to prove that the total variation is uniformly 
bounded by 0{l)5i. As we will see, the crucial point is to prove that, if Tot Var{u2:(cr)} < 0(1)^1 
for all a < t, then it holds an estimate of order two on the integrals of the source term: 

/ [ \Si{cr, x)\dxdcr < 0{l)Sl. (1.24) 
JO Jo 

To show ((1.2411 we will basically deal with each of the term that appear in the expression of si 
separately. Some of the estimates are based on the same techniques described in 0: in particular 
we will use the interaction, area and length functional introduced in the boundary free case. Some 
estimates, on the other hand, require quite long computations and can be found in the appendix. 
In Section [S] we will prove the stability of the vanishing viscosity approximation with respect to 
perturbations. More precisely, let mq, Ubo, u^l and vq, vto, VbL be the initial and boundary 
data of two solutions u and v of problem 1(1.16(1 : we will show that there exists a constant Li such 
that 

\\u{t) - w(t)|lLi(o,L) < il(||wo - Wo||l1(o,L) + ||wbo - Wbo||Li(0,i) + ll^bL - Wf,L || Li (0, t)) ■ 

Moreover, one has also stability with respect to time: if u is a solution to ((1.16(1 then 

\\u(t) - ii(s)||ii < L2{\t - s| + |Vt - Vi|) 
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for a suitable constant L2. We will see that the constants Li and L2 depend uniquely on the matrix 
A and on the bound 5i on the total variation of the initial and boundary data. We will actually give 
just a sketch of the proof of the stability, since we will show that one can employ the same tools 
used to prove the BV estimates and repeat with minor changes the computations of Sectional 

One can then get back to the solution u"^ of the original problem and obtain that for all 

e > it satisfies 

TotVar{u^(t)} < ©(l)^! Vt > \\u'' {t) - u*\\oo < 0{l)5i yt>0 

\\u^{t) ~ w''(t)||Li(0,L) < Li(\\uo - Wo||l1(0,L) + \\ubO - W6o||Li(0,t) + ll'^bO - "bL || Li (0, t)) (1-25) 

\\u%t) - u%s)\\li < L2{\t ~s\ + V^\Vi-^s\). 

In the last estimate, uq, Ubo utL and vq, Vbo vtL are the initial and boundary data for two solutions 
and of (fO|) . 

The uniform bound on the total variation of the solutions of 1)1.3(1 ensures that for any 
(mo, Ubo, Ubi) € 1^0 X i^b X Zifc, for any t > and £„ — > 0+ there is a subsequence e„^. such that 
(t) converges in L^{0, I) to some limit function we will denote by pt[uQ, Ubo, Ubi]- Letting e — > 0"*" 
in H1.25|l one finds that the limit satisfies the stability estimate 

Pt[uo, UbO, Ubl] -Ps[vo, VbO, Vbl] ^ <Lii \\vo - Uo\\l^{o,1) + \\vbo - "6o||l1(0,+oo) 

^ (1.26) 

+ \\vbl - Ubl\\Li{0,+oo) ] + L2\t - s\. 

By a standard diagonalization procedure one can show that there is a subsequence that converges 
for any rational time t and for any {uq, Ubo, Ubi) in a countable dense set oUAq xUb y-Ub', the density 
is here intended in the norm. Then by the estimate H1.26|l pt[uo, Ubo, Ubi] must be defined on 
close sets of times and boundary and initial data. Hence [mq, UbQ, Ubi] is defined for any t > and 
for all (-So, "Sfco, Ubl) eUa xUb x Kb- 

One can actually check that the operator 

S : [0, +oo] X Uq X Kb >^ Ub ^ Vq X Kb X Kb 

{t, Uo, Ubo, Ubl) PtN, Ubo, Ubl], Ubo{ - +t), Ubi{ - +t) 



satisfies the semigroup property 

To complete the proof of Theorem ll.2l one is therefore left to show the uniqueness of the semigroup 
of vanishing viscosity solutions: indeed, different sequences u'^"(t) could a priori converge to 

different limits. 

The proof of the uniqueness of the vanishing viscosity limit can be found in Section 16.41 and, 
following the same ideas as in Q, the crucial step will be to show that the semigroup defined via 
vanishing viscosity approximation is actually a viscosity solution in the sense of 0. 

We refer to Section for the precise statement, here however we underline that the definition 
of viscosity solution is based on local estimates that ensure, roughly speaking, a "good behavior" in 
comparison with the solutions of a suitable Riemann problem and of a suitable linear problem. 

The notion of viscosity solution was first described in the conservative boundary free case in 
|1()) and was strictly connected to the definition of Standard Riemann Semigroup (SRS) that was 
introduced in the same paper. For completeness, we recall here that a SRS is Lipschitz continuous 
with respect to the norm and in the case of piecewise constant initial data locally coincides with 
the standard Riemann solver defined by Lax in |2B]. In TU' it is proved that if a SRS semigroup 
exists, then it necessarily coincides with the wave-front tracking limit and with the viscosity solution. 
One of the main advantages one gains introducing the notion of viscosity solution is therefore the 
characterization of global behaviors through local ones. 

The definition of SRS semigroup and of viscosity solution was extended to conservative boundary 
value problems in 0. Moreover, in the same paper it was proved that, also in the case of an initial- 
boundary value problem, if a SRS exist then it necessarily coincides with the wave-front tracking 
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limit and with the viscosity solution. Hence the uniqueness of the SRS semigroup comes from the 
uniqueness of the wave- front tracking limit, proved in |20| . 

From the previous works it is clear that a crucial step in the definition of viscosity solution is 
the description of the Riemann solver and of the boundary Riemann solver. 

As mentioned before, a solution of the Riemann problem in the boundary free case was introduced 
by Lax ([2b\) for conservative systems in the case of linearly degenerate or genuinely non linear 
fields. Such a definition was then extended by Liu ( 27 ) to very general conservative systems. The 
characterization of the Riemann solver for non conservative systems was introduced in jTj, where it 
was also proved the effective convergence of the vanishing viscosity solutions and it was extended in 
the natural way the notion of SRS and of viscosity solution. 

As concerns boundary Riemann solvers, a solution of the initial boundary value problem 

ut + A{u)ux = 

(1.27) 

u{t, 0) = Ub u(0, x) = uo, 

was proposed in |21j in the case of systems in conservation form with only linearly degenerate 
or genuinely non linear fields: such a boundary Riemann solver is in general different from the 
one defined by the vanishing viscosity limit (some more precise considerations can be found in 
Remark 16.111 . On the other side, in [251 131L IT] and |2| it was considered a quite general boundary 
condition, which turns out to be compatible with the one defined by the limit of vanishing viscosity 
approximations: we refer again to Remark lfi.ll for a more precise statement. We underline, moreover, 
that a study of the boundary conditions defined by the limit of the general parabolic approximation 

can be found in |28l 1221 1281 1291 IHH) in the case of systems in conservation form. Finally, the 
Riemann solver for boundary value problems non necessarily in conservation form was first described 
in P); in this paper it was also extended in the natural way the notion of SRS and of viscosity solution. 

In Section 16.11 we will describe the Riemann solver and the boundary Riemann solver defined 
by the vanishing viscosity limit, which however have an interest in their own. The problem dealt 
with is actually a particular case of the one solved in 0, where also the characteristic case was 
considered, but since the reduction to our case is not completely trivial, we will describe it explicitly. 
In particular, we will consider the vanishing viscosity solution of the boundary Riemann problem 
(|1.27|) . Let ^(0+) = lim2.^Q+ u{t, x) be the trace of the solution on the axis a; = 0, which does not 
depend on time since the solution u is self-similar. We will show that there exists a solution of the 
ODE 

A{U)U, = U,, (1.28) 

such that 

U{0) = Ub, lim U{x)^u{0+). 

X — >-\-oo 

In other words, the boundary datum Ub does not necessarily coincide with the trace u(0"'"), but it 
certainly lays on the stable manifold of u{0~^) with respect to the ODE (|1.28|) . 

Remark 1.3. The fact that the bounds on the total variation are uniform with respect to the length 
L of the interval implies that, for any fixed e > 0, one can let L +oo in (|1.16(l . Hence, coming 
back to the original system (|1.3I) one finds that also the solutions of 

wf -I- yl(u^)Mj = eu|j., a;e]0, -foo[, te]0,+oo[ 
< u'^{0,x) = uo{x) 

^ u%t,0)^uba{t) 

have total variation uniformly bounded with respect to e. 

Hence the analysis of the vanishing viscosity approximations of the initial-one-boundary value 
problem can be deduced as a limit case from the study of the two boundaries case. 
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2 Parabolic estimates 



In this section we will find a representation formula for the solution to (|1.16fl 

ut + A{u)ux ^ Uxx, x g]0, L[, i e]0, +oo[ 

u{0,x)^uo{x), (2.1) 

u{t, 0) ^ Uboit), u{t, L) ^ UbLit) 

with initial and boundary data satisfying (|1.7|) . I|1.8() and l|1.9|l . The aim is to prove that the solution 
of (|2.1|l is regular and that the norm of the second derivative ||ua:a;(t)||ii(o.L) is bounded, as soon 
as the total variation of u{t) remains small. We will regard (|2.1|l as a perturbation of the linear 
parabolic system with constant coefficients 

ut + A*Ux - Uxx = 0- (2.2) 
Here and in the following we will assume A* = A{u*) and A* = Xi{u*). 

2.1 The convolution kernels 

The fundamental step is to study the equation (|4.1|l in the scalar case, because the Green kernel for 
the general vector case ()2.2|l follows by using the base of eigenvectors of A* . Thanks to the linearity, 
we split the Green kernel of the equation 

zt + X*Zx - Zxx = (2.3) 

into 3 parts: 

1. A^i (t,x,y) is the solution of (|2.3|) with zero boundary conditions and initial condition 

A^HO, x,y) = dy z/e]0,L[. 

This function is given by 

A^* (t, y) = ( G{t,x + 2mL-y)-G{t,x + 2mL + y)](j)^'{t,x,y), (2.4) 



m — —oo 



where G{t, x) — (e ^^/^*)/2>/7rt is the standard heat kernel and 



~^'it,x,y)=cxp(^ix-y)-^t 



2. J^^"{t,x) is the solution of with zero initial datum and boundary conditions 

J^'^°it, 0) = 1 J^'°(t, L) = 0. (2.5) 

It follows that 

J^^°{t,x)=Aexp{X*x)+B-J^ A^'{t,x,y)(^Aexp{X*y)+B^dy, (2.6) 

with 

A = ^ B=^^ . 
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3. ^{t,x) is the solution of H2.3|l with zero Cauchy datum and boundary conditions 

J^'^(t,0) = J^'^(t,L)^l (2.7) 

and it is given by 

J^'^^{t,x)^Cexp{X*x)+D- A^*{t,x,y)(^Cexp{X*y) + D^dy, (2.8) 

where 



Note that aU the coefhcients A, B, C, D remain bounded as L ^ +oo. Moreover, one can ap- 
ply the maximum principle and, via a comparison with the constant solutions, finds that < 
J^' °(t,x), J^' ^{t,x) < 1. Hence the integrals 

J^'°{t,x)v'{t)dt [ J^'°{t,x)v'{t)dt 



are well defined for every function v(t) € BV{0, +oo) and for every T. 

In the following, we will also need a further convolution kernel A^* {t, x, y) such that 

Ky^{t, X, 2;) + Af (t, X, y) = 0, 

i.e. 

A^'it, x,y)^ / A^' {t, X, z)dz. (2.9) 

To get the previous formula we have arbitrarily imposed A^» (t, x, L) = 0. 

Note that A'*'^ (t, x, 0) is the derivative with respect to a; of a function z(i, x) which satisfies 

z{t, x) + J-^' °(t, x) + J^' ^{t, x) = 1. 

Hence, 

A^- [t, X, 0) + J^* °(<, x) + J^' ^(t, x) = 0. (2.10) 

The following proposition provides some basic estimates on the convolution kernels we will need 
later. 

Proposition 2.1. The convolution kernel A"^' satisfies 

||A^*(i, y)|Ui <0(1) ||A^'(t, y)|Ui <0(1)/Vi Vt< 1, ye]0, L[. (2.11) 
The following estimates hold for the boundary kernels J^' ^ , J^' ^ : 

< J-^- °(t,a;), J^^^(t,x) < 1 yt>0, X e]0,L[ 

\\J^'\t)\\L^,\\J^'\t)\\L^ <Oa) VO<t<l, (2.12) 

I|./..1°(0IIli,I|JxI^(0IIli <o(i)/\/i V o<<< 1. 

T/ie auxiliary convolution kernel A^' satisfies estimates analogous to those of A^> : 

||A^'(i, y)|Ui <0(1) ||A^*(t, y)|Ui <0(1)M VO<t<l, ye]0, L[. (2.13) 
The proof of the proposition can be found in the Appendix I A. 1 . ll 

Now we are ready to deal with the vector case. Let r* , I* i — 1,2 be respectively the left and 
the right eigenvectors of A* = A{u*). We define the matrix kernels 

A* := ELi A^*r* ® /*, A* :^ ^Li ^'^r* ® Z*, 

(2.14) 

J*o ,^J2i=iJ^^"r*(g>l*, J * ^ Ell ® • 

By construction these are the matrix kernels for the initial data corresponding to the cases 1 , 2 and 
3 considered above (equations 12. 4|) . 1)2.5(1 and ((2.7|l respectively). 
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2.2 Parabolic estimates 

The solution of equation (|2.1|l can be written as 

u{t,x)= I A*{t,x,y)uo{y)dy + uo{0)J*"{t, x)+ [ J*"{t - s,x)u'^Q{s)ds + uo{L)J*^{t, x) 
Jq Jo 

+ I J*^{t - s,x)u',,i^{s)ds + / / A*{t - s,x,y){A* ~ A{u))uy{s,y)dyds, 
Jo Jq Jo 

(2.15) 

and therefore, recaUing H2.10|l and integrating by parts, 



ix{t, x) = I A*{t,x,y)uQ{y)dy + J*°{t - s, x)u',,„{s)ds + J* {t - s , x)u'f,L{s)ds 

Jo Jo 

A*{t-s, X, y)(^(^A* - A{u)Juyy - DA{u)[uy (E) uyj^ {s, y)dyds (2.16) 



Jo 

+ j;^{t, x){uo{L)~uo{0)) - / + J;^)(i-s, x){A* - A{u))u^{s, 0)ds. 

Jo 

From the previous expression we immediately have that, as long as it can be prolonged, the solution is 
regular. Moreover, the local existence of a solution of equation H2.1|l follows from the representation 
formulas (|2.15|l and H2.16|l via the contraction map theorem. 

We can now use the representation (|2.16|) to prove the following proposition. 

Proposition 2.2. // \\u^{t)\\Li < 0{l)Si for all t e [0, 1], then 

\KAmL^<^^ Vie [0,1]. 

Proof. From (|2.16|) we get 

pL nt pt 

Uxx{t,x)^ I Al{t, X, y)u'o{y)dy + J*°{t - s,x)u'^Q{s)ds + J^J^ {t - s, x)u'^^{s)ds 
Jo Jo Jo 

J A* (i - s, X, ?/) ^ (^A* - A(u)^ Uyy - DA{u) (^Uy (g) ^ (s, y)dyds 

+ j;/(t, x){u^{L) - 7.0(0)) - 1* {J:° + j;/)(t - s, x){a* - A{u))u.,{s,Q)ds. 

(2.17) 

The previous representation formula shows that the function 1 1— > |lwxa:(i)|lLi is continuous. 
We claim that there is a constant C independent from L such that 

Indeed, for a fixed large constant C, define 

T = inf|i: \\u,j,^{t)\\ > -jzSi 

The time r is strictly bigger than if C is sufficiently large, since by hypothesis ||moI|z,i is finite. 
Moreover, one has ||wxa;(''')||Li = C5i/\/t thanks to the continuity of the map t ^ \\uxx{t)\\L'^- 
From ifTTTjl it follows that 

\\uUr)\\L^ - < II A;(r)|Ui ||u^,||l^ + 0(l)<5i T \\uyy{s)hi \\A*,{t - s)h^ds + 25, f ^EL ds 

Jo Jo VT - s 

0(1) X , r 0(1) c .rf. 



p. ^251 r 

Jo 



Vs{t- s) 
< + 20(1) CSl + 20(1) V^Ji, 
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which is a contradiction if C is large enough and 5i sufficiently small. In the previous estimate we 
have used the bounds 

IkUlk- < < Si r J- . ds = vr. 

Jo y's[T-s) 

□ 

If t > 1 and ||Ma;(s)|lLi < 0{l)Si for any s G [0, t] , we can apply the previous proposition to the 
interval — 1, i] and obtain 

\\u,x{t)\\L^ < 0{l)Si t>l. 

Since the derivative Ux is regular, this implies in particular that, if ||u2;(s)||ii < 0{l)di for any s <t, 
then ||wa;(t)||L°° < 0{l)Si if t > 1: in other words, as long as Ux remains small in the norm, it 
remains small in the L°° norm too. 

3 Gradient decomposition 

3.1 Double boundary layers and travelling waves 

In this section we will introduce a suitable decomposition of the gradient of the solution to H1.16|l . 

ut + A{u)ux = Uxx, X e]0, L[, <e]0, +oo[ 

< u{0, x) = Uq{x), 

u{t,0) ^ Ubo{t), u{t, L) ^ UbL{t). 
We will employ a decomposition in the form 

Ux ^ Vifi + V2f2 + Plfl + P2r2, (3.1) 

where the first two terms correspond to derivatives of travelling waves and the last two correspond 
to the derivative of a double boundary profile. More precisely, pi is the part of the double boundary 
profile exponentially decaying as a; — > +oo, p2 is the part exponentially decaying as a; — > — oo. 

The principal results of this section are the construction of the vectors fi, f2, the description of 
a decomposition of Ux in the form 1)3. l|l . the computations of the equations for the 4 components 
t^i, V2, Pi, P2 and finally the choice of the boundary conditions for the same components. In the 
description of the decomposition we will focus mainly on the construction of the double boundary 
profiles, because the construction of the travelling wave profiles follows the same steps as in |7j. 

The construction of the double boundary profile is based on the following idea: in the linear case, 
one finds that there is a solution of the boundary value problem 

{Ux 
Px = A{u)p, (3.2) 
u{0) = Ubo, u{L) = UbL 

and such a solution is the sum two components: one exponentially decaying as a; ^ +oo, the other 
as a; — + — oo. Moreover, when the length L is very large the solution has the behavior illustrated 
in figure 121 (on the left): it is very steep near the boundary a: = because of the presence of the 
exponentially decreasing component, then it is almost horizontal in a large interval and then it is 
steep again near the boundary x — L because of the presence of the exponential decreasing part. 

The idea is to try to simulate such a spatial behavior also in the non linear case: in this way, when 
L is large enough the derivative of the double boundary profile is concentrated near the boundaries 
a; = and x = L and therefore there is essentially no interaction with the travelling wave profiles 
inside the domain. This behavior is the same one observes in the hyperbolic limit, where in ]0, L[ the 
solution is generated only by travelling wave profiles. We will find out that, if |C/bo ~ ^-^6l| is small 
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enough, then there exists indeed a solution of the boundary value problem (|3.2I) with the behavior 
illustrated in figure |21 

In this way, we construct the functions pifi{u, pi) and P2'^2(w, P2)- however, since the decom- 
position H3.1|l is a 2-dimensional vector equation in 4 scalar unknowns, we have some freedom in 
assigning the initial and boundary data for vi, V2, Pi and p2. The detailed description of the bound- 
ary conditions can be found in Section 13.31 but the crucial idea is to impose some conditions that 
allow the component pi and p2 to behave like the derivative of a double boundary layer, and thus 
to be independent from the choice of the initial datum and to be concentrated near the boundary 
X = OT X = L, respectively. On the other hand, we want to impose some conditions on the compo- 
nents vi and V2 that forces them to behave like the derivative of waves in the hyperbolic limit, thus 
flowing out from the domain through the boundary x = (waves of the first family) or through the 
boundary x = L (waves of the second family). 

Moreover, we have also some freedom in assigning the source terms, as it will be clear in Section 
13.21 again the basic idea we will follow is that p2 , which corresponds to the component of the double 
boundary profile exponential decaying as x —^ +oo, should be affected only by the datum in x = L. 
Since in general the source term are spread on the whole interval ]0, L[, we will impose that the 
equation for p2 has no source term. 

3.1.1 Double boundary profiles 

As a first step, we characterize the solutions of the system 

1"^=^,,, (3.3) 
{ Px = A[u)p 

that converge with exponential decay to some value (u, 0) with u in a small enough neighborhood 
of the value u* defined by the relation (|1.8() . Since {u*,0) is an equilibrium point, we can consider 
the linearized system, whose center and stable subspaces are given by 

V'^^Ip^ 0}, V = span(ri(w*)), V = span(r2(u*)). 

Let {pi, P2) be the coordinates of p with respect to the base defined by the eigenvalues ri{u*) 
and r2{u*) of A{u*): thanks to the center-stable manifold theorem, there exists a regular function 

4):{{u,pi) ■.\u~u*\,\pi\<e} QV" (BV R 

{u,pi) K-* p2 = (f>{u,pi), 

which parameterizes the solutions of (|3.3() that do not blow up exponentially for x +00. In 
our case, one can see that this manifold is made by the orbits which converge for x — s- +00 to an 
equilibrium (u, 0), with u close to u* (figure^. In particular this manifold is unique. 

The dimension of this manifold is dimV^ + dimV^, i.e. 3 in our case. Since pi — implies 
P2 — 4'{u^pi) = 0, we can set (j){u^pi) ~ pih{u,pi) and M'^'^ can be described by the following 
condition: 

p^piri{u*)+pih{u,pi)r2{u*) =^pi(^ j.^ ^ := pifi(w,pi). 
Inserting the previous expression in the system (|3.3I) . one obtains 

A{u)pifi = (pi^i)^ =Pixri + (pi)^Dfifi +pipixfip. 
Let £1 — (1, 0): if we multiply the previous expression by £1 we obtain, since A is triangular, 

AiPi = Pix, 

and hence 

A{u)piri = Aipifi + (pi)^Dfifi -I- Xipfrip. (3.4) 
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Figure 1: the center-stable manifold A^'^* and the center-unstable manifold with orbits expo- 
nentially decaying to an equilibrium point as a; — > +00 or a; —00, respectively 




It follows that 
and therefore 



fi(ii, 0) = ri{u) V u, 



In a similar way one can also define a regular, 3-diniensional center-unstable manifold Al"^" 
containing all the orbits that as a; — > —00 converge with exponential decay to some point (u, 0) with 
u close to u* . The manifold is parameterized by V^® V"; moreover, since the matrix A is triangular, 
one can choose 

' 



h = r2{u) = 



1 



The manifold Al^" is thus described by the relation p = P2r2- 

As a second step, we show that the functions pivi and P2''2 indeed allow us to construct a solution 
of the two-boundaries value problem 



Decomposing Zx as 

Zx = Pifiiz, pi) +P2r2 

and using the relation (|3.4|l . we obtain the system 

Zx = Pifiiz, pi) +P2r2, 

Pix = Ai(z)pi, 

P2x = A2(z, Pl)p2 



(3.5) 



where we have defined 



A2(m, Pi) :== >^2{u) Dfir2), 



(3.6) 



(3.7) 
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where the vector ^2 satisfies {£2, fi) =0 and {£2, ^2) = 1. Hence, while in the hnear case the two 
components of the solution of the system H3.5|l are decoupled, in the general case there is a coupling 
in the equation of z, and in the choice of A2, which is in some sense the effective eigenvalue for p2. 
Note that 

\X2iu,pi)-X2iu)\<0{l)pi. (3.8) 

An application of contraction principle ensures that, if|C/f,o^C^f)L| l£ ^1 for a small enough 5i, then 
the above system with boundary data z(0) = Ubo, z{L) = UbL has a unique solution. Moreover, 
one also finds that |A2(u, pi) — \{u)\ < 0(1)^1. 

Since Ai < 0, A2 > for di 1, we obtain that pi is exponentially decaying, while p2 is 
exponentially increasing. We can thus figure the double boundary profile as follows (figure EJ: when 
the length L of the interval is very large, the solution will be steep near zero, because in that region 
Pi varies exponentially fast. Then it will be almost horizontal for a long interval and becomes again 
very steep in a left neighborhood of x — because p2 increases exponentially. 



3.1.2 Travelling waves 

We refer to .7 for an exhaustive account of the analysis that allow the definition of the decomposition 
along travelling waves: here we will only recall for completeness the crucial steps. 
Consider the system 

{Ux ^P 
p, = {A{u) - al)p (3.9) 
a, =0 

and an equilibrium point {u* , 0, \i(u*)). The center manifold theorem ensures that the center space 
= {p = 0} parameterizes a center manifold M.'^. This manifold contains all the solutions of H3.9|l 

that do not diverge exponentially neither as x — > — cxd nor as a; ^ +00. 

It can be shown that the center manifold Al"^ around the equilibrium (m*, 0, \i{u*)) is described 

by a function Vifi(u, Vi, ai). Since A is triangular, one can take 

ri(u, V,, ai) = ^(^^ ) , f2(«, V2, ^2) ^ ( J ) , 

for some suitable function m (in general different from the function / in the vector fi). One can 
moreover show that the following equations hold: 

A{u)fi = Aifi + wiDrifi + wi(Ai — cri)fii,, 
fi(u, 0, CTi) = ri(M) Vu, ai, |r(u, ui, cti) - ri(w)| = 0(l)wi, fi„^0{l)vi. 

Here and in the following we will denote by (^i, £2) the dual base of (fi, r2). 
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3.1.3 Gradient decomposition 

We set 

Ux ^ vihiu, vi, ai) + V2r2 +Pifiiu,pi) +P2r2 / \ 

ai = Xi{u*)-9i^ + Xi{u*)y (3.10) 
ut = wiri{u, vi, ai) + W2r2 ' 

The function Q is here and in the following an odd cutoff such that 

r s if|si <<5 

Q{s) = l if|s|>3(5 5i«<5<g. (3.11) 

I. smooth connection if (5 < s < 3(5 

The choice of the speed a follows from the analysis of the boundary free case, [71. 

Note that (|3.10|) is a system of 4 equations in 6 unknowns: as we underlined in the introduction, 
this will allow some freedom in choosing the boundary conditions for u;, « = 1,2 and pi, i = 1,2. 
More precisely, we will proceed as follows. 

1. We will insert H3.1()|l in the parabolic equation (|1.16(l . This will generate a system of 4 equations 
in 6 unknown. 

2. We will obtain the equations for u^, wi, pi, i = 1, 2 by assigning in a suitable way the terms 
obtained. 

3. We will impose boundary and initial conditions on each of the 6 equations obtained. This 
procedure selects one and only one solution for each of those equations. 

The decomposition (|3.1U|I is thus complete. We observe that the idea is to let the equations to choose 
the components in the decomposition, by only imposing reasonable initial-boundary conditions and 
by assigning carefully the terms obtained by inserting (|3.1()(l in the system (I1.16|l . 

3.2 The equations satisfied by v^, p^, Wi i = 1, 2 

These equations are obtained via the computations in Appendix IA.2.11 inserting the components 
Vi, Pi Wi, i ^ \, 2 \n the equation 

Ut + A{u)Ux - Uxx = 0, 

we find 

Vit + (Ait;i):r - Vixx + Pit + {XlPl)x - Plxx = 

V2t + {X2V2)x - V2xx + P2t + {X2P2)x ~ P2xx = Sl(i, x) 

Wit + iXlWi)x - Wixx = 

W2t + iX2W2)x - W2xx = hit, x) 

for some function Si{t, x) i = 1, 2 whose explicit expression can be found in the appendix. Moreover, 
as it is shown in the Appendix IA.2.11 from the equation 

Ut = Uxx - A{u)ux 

one gets the relations 

Wl = Vix - XlVi + pix - XlPl 

(3.12) 

W2 = V2x - X2V2 + P2x - X2P2 + e{t, x) 

for a suitable error term e{t, x). The following Proposition (whose proof can be found in the 
Appendix I A. 2. 2|l gives the form of the source terms: 
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Proposition 3.1. The following estimate holds: 

a:)|, \~S2(t, x)l |e(i, a:)| < 0{l)\ ^ [\v,\(\v,\ + + + \w,.,\) + \w,\(\w,\ + ^ 

l(|P + IP2I 

2 

X{|t«i|<<5i|t.i|} + \wi +(^ivi\{\vi\ + \vi.x\ + \wi\ + 



(3.13) 



Following the denomination of P], we will denote the above terms as follows: 
1. interaction between waves of family 1 and family 2 



Vj \ + \Vjx \ + \Wj \ + \Wjx\ ] + 



\w^\{\wj \ + \v^ 



3x1 h 



2. interaction of travelling waves with boundary profiles 

Vj \ + \Vja:\ + \Wj\ + \WjM; 

hi 

3. interaction among boundary profiles 

\Plx - AiPi|(^|p1:e| + \P2 

4. CTi is not constant 



\wivix - viwixl + vl 



X{iu.ii<5ii«il}; 



5. the cutoff fimction 9 is active 

\wi + aivi\(jvi\ + \vix\ + \wi\ + \wix\j. 

Since the component p2 of the boundary profile should remain close to the boundary x — L, and the 
source si is in general spread in the whole interval [0, L], we split the previous expression as follows: 



Vlt + {XlVl)x - Vixx = 

V2t + i\2V2)x - V2xx = Si{t, x) 



Pit + i>^lPl)x ~Plxx = 
P2t + i>^2P2)x ~P2xx = 



3.3 Boundary conditions 



To conclude the characterization of the equations satisfied by Vi, pi, Wi, we have to assign the 
boundary conditions. The basic idea is that each component Vi, pi, i = 1, should behave like 
a travelling wave or a boundary profile, respectively. More precisely, we can make the following 
observations: 

1) In order to behave like a double boundary profile, pi and p2 should be independent from the 
initial datum, hence we are led to impose 

pi{0,x) = 0, p2(0, x) = 0. 

It follows that the initial data for vi and V2 are given by 

vi{0, x) = {h, u'oix)) V2{0, x) = {£2, u'oix)). 
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2) To emulate the behavior observed in the hyperbolic limit, the waves of the first family should 
disappear when hitting the boundary a; = 0, and the waves of the second family should disappear at 
X — L. To understand what kind of boundary condition it is convenient to impose, one can observe 
that an integration by parts leads to 



— J \vi{t, x)\dx =^ J signal - Aii^i^ dx 



Sv=o{vix - Xivi)dx + 



sigm;i(wi^ - XiVi) 



< 



signal (dij: - Xivi) 



f \v2{t, x)\dx < f f \si{s, x)\dsdx 
di Jo Jo Jo 



Signw2(w22; - A2U2) 



(we have used the inequality S.^^o'^^x < 0). To minimize the increment of ||wi(t)|lLi due to the 
interactions with the boundary we impose 

0) = 0, V2{t,L) = 0, 

and integrating with respect to t the previous equations we get 



vi{t, x)\dx < / \vi{0, x)\dx + / \vix ^ XiVi\{s, L)ds 



\v2{t, x)\dx < / \v2iO, x)\dx + 
Jo 

We have used the following observations: 

wi(0) = = 
V2{L) = . 



(3.14) 



\si{s, x)\dsdx + / \v2x - X2V2\is, 0)ds. 



lim sign(t;i)wia;(a;) > 

x^0+ 

lim_ sign(t;2)i;2a;(a;) < 0. 

X — »L 



(3.15) 



If one inserts the previous Dirichlet condition on vi i = 1, 2 in the decomposition H3.1()|l . obtains 
the foUowings boundary conditions for pi : 



Pi{t, 0) = {£1, ux{t, 0)), P2{t, L) = {i2, ux{t, L)) - pi{i2, h). 



(3.16) 



3) Since pi should be located near x = 0, and p2 near x — L, we would like to impose that the 
increment of due to the datum at a; = X is minimal, and similarly that the increment of 

||P2||li caused by the boundary datum in x = is as low as possible. Since the values pi(t, 0) and 
P2it, L) are already determined, we will impose on pi some condition at x = L and on p2 at x ^ 0. 
We observe that an integration by parts like the ones performed before leads to 

/ \pi{t,x)\< / \pix- Xipi\{s,0)ds+ j \pix - Xipi\{s, L)ds. 
Jo Jo Jo 

Hence we are led by the previous considerations to impose 

[pix - Xipi){t, i) = 0. 

Similarly, we impose 



{P2x - X2P2){t, 0)=0. 

From these two equations we obtain the boundary conditions for i>i, V2'. indeed, we have 



(3.17) 
(3.18) 



and 



(wi,x-Ait;i)(t, L) = (£1, ut{t, L)) 
(V2,x - X2V2){t, 0) = (£2, ut{t, 0)) - e{t, 0). 



At this point, the initial-boundary data are perfectly determined for all the components Vi, pi 
i = 1, 2, and thus the decomposition is complete. 
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4 BV estimates 



Aim of this section is to prove the fohowing theorem, which constitutes the first part of Theorem 

o 

Theorem 4.1. Let u{t, x) be the local in time solution of the 2x2 system 

Ut + A{u)ux = U^x 

u(0, x) = uo(a;) (4.1) 

u(t, 0) = Uho(i) u{t, L) = UbLit) 
and suppose that the boundary and initial conditions are regular and satisfy 

< Si k = 1, . . .n, 



d'^uo 




d^Ubo 




d'^UbL 




dx^ 




dt^ 


Li(0,+oo) 


dt^ 


Li(0, +oo) 



for some 5i sufficiently small. 

Then u{t, x) is defined V t > and its total variation is uniformly bounded: 

\\Ux{t)\\L^O,L) < CSi 

for some constant C which does not depend on L. 



(4.2) 



It is enough to prove that there is a constant Sq such that kSi < So << 1 with k small enough 
and such that the following holds: if di is small enough and ||'«a;(s)|jii < C6i Vs G [0, t] then 



t rL 



\Si{<J, x)\dxda- < 0{1)S^, 



O^JO 



t nL 



Ihia, x)\dxda- < 0{l)Sl, 



O^JO 



\v2x - >^2V2\icr, 0)da < mSi, / \vix - XiVi\{a; L)d(7 < mSi, 
\pix - >^iPi\{(J, 0)d(T < mSi, / |p2x - ^21(0", 0)da- < TO(5i, 



(4.3) 



for some constant m that does not depend on C. 

Indeed, suppose the previous implication holds. From the representation formula H2.16|l it imme- 
diately follows that the function 1 1— > ||wa;(t)||^i is continuous: hence, it will satisfy ||ua;(t)||^i < CSi 
if t is small enough, since the total variation of the initial datum is bounded by 5i. 

Suppose by contradiction that r is the first time such that ||u2:(T)||ii = CSi. Then we use the 
equations 

Vit + iXlVi)x - Vixx =0 pit + {XlPl)x - Plxx = 

V2t + iX2V2)x - V2xx = Si{t, x) P2t + {X2P2)x - P2xx = 

and the boundary conditions described in Section 13.31 and, integrating by parts, we get 



/ \ux{t, x)\dx \v,{t,x)\+ \p,{t, x)\dx \v^{0,x)\ + 

Jo „_i Jo Jo „_i Jo Jo Jo 



\si{a, x)\dxda- 



+ \v2x - X2V2\ {a; 0) da + \vix - XiVi\{a, L)da + \pix - Xipi\{a, 0)da 
Jo Jo Jo 



+ / \P2X - X2P2\{(T, L)da < {Am + 2)Si + 0(1)31 < CSi 



if C is large enough: this contradicts the assumption ||u2:(T)||ii = CSi. 
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Note that since all the functions in the right hand side of H3.13|l are continuous (and hence 
bounded on [0, L]), we have that 

L 

\s^{a, x)\dxda < 0{l)6i i = 1, 2, (4.4) 



Jo 

for s small enough. Hence to prove 1)4.3(1 we can suppose that 1)4. 4|l holds for any s G [0, t]: since we 
will show that actually 

/ [ |Sj(cr, x)\dxda < 0{l)6l, i ^ 1, 2, 
Jo Jo 

the assumption will be a posteriori justified since k6i < 6q << 1. 

We will proceed as follows: in Section 14.11 we will show some elementary estimates, while in 
Section ^21 we will introduce suitable functionals that allow the estimates 

J2 (I^^KI^il + I^J^I + l^jl + I^'J^I) + l^'^KI^il + \vjx\)){crx)dcrdx < 0(1)(5?, 



Jo 



t i-L 

\wivix - ViWix\{cr, x)da < 0{l)dl, 



Jo 

t 



Jo 



^2 1 Wl 

Vl 



2 



X{\wi\<5i\vi\}i(^, x)dadx < 0{1)S 



In Section Ol we will consider the term 

rt i-L 



\wi + (Tiwi|(|wi| + \vix\ + \wi\){a, x)dadx, 



and prove a bound of order Sf 



JQ 
2 



4.1 Elementary estimates 

This section is devoted to the estimates which can be obtained by elementary techniques, like the 
maximum principle. We will in particular show that the components pi, i = 1,2 are exponentially 
decaying as one moves far away from the boundary, and that their decay exponent does not depend 
on the interval length L. Moreover, by introducing various functional, we estimate the boundary 
data assigned to the components vi, V2 and prove that the functions Vi are integrable along all 
vertical lines {x = const}. This means that, as in the boundary free case, the profiles of travelling 
waves just cross the vertical lines. 

4.1.1 Estimates via maiximum principle 

We will first deal with pi. The results in Section ensures that 

\\u,{t)\\L^ < \\u,,{t)\\Li < 0{l)Si. 

Hence it follows that 



\pi{t, Q)\ = \{h, ux{t, 0))| < k6 



for some k large enough. 

The equation satisfied by pi is 

Pit + Ai(u)pi3; + Xix{u)pi - pixx = 0. 

This is a linear equation, with coefficients depending on the solution u(t, x). Let 2c be the separation 
speed defined in (|1.5|) and 

q{x) — kSi exp ( — ca;/2). 
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Since |Aia;| < 0{l)Si and 5i « 1, q satisfies 

qt + Mqx + Mxq - qxx > 0. 

Hence the difference {q — pi) satisfies 

{q-pi)t + M{q-pi)x + Mx{.q-pi) - {qxx -pixx) > 
{q~Pi){t. 0)>0 
(q - Pi) - Ai(q - pi)^ ] {t, L) > 0. 



By standard techniques it foUows that (q — pi)(t, x) > for any t, x and hence 

\pi{t, x)\ < kSi exp(-cx/2). (4.5) 
The boundary condition on p2 satisfies the following bound: 

\P2it, L)\ = \{i2, u,{t, L)) -plih, h)\ < o{i)Si, yt, X. 

Since \pi{t, x)\ < kSi, then from H3.8|l it follows that IA2 — A2I < 0{l)Si and hence in the same way 
as before one can prove 

\p2{t, x)\ < 0{l)Si exp(c(a; - L)/2), Vi, x. (4.6) 

From (|4.5|l it follows 

\\pimL^<0{l)Si, \\viit)\\Li <0{l)Si 
and, since < 0{l)6i, 

llfilli- <0{l)Si. 

Analogously, from H4.6(l it follow 

lb2(t)IUi < o(i)<5i, IMtjU^ < o{i)Si, \\v2\\l^ < o(i)(Si. 

The following proposition summarizes the results obtained in this paragraph: 

Proposition 4.1. Let pi, Vi be the solutions of 11.21|l with the boundary conditions described in 
Section\^ Then 

\piit, x)\ < 0{l)Si exp(-c.i;/2), \p2{t, x)\ < 0{l)Si exp(c(a: - L)/2), 

where 2c is the separation speed defined by (|1.5|l . 
The previous estimates imply 

\\Prit)\\L^ <0{l)Si, \\v,it)\\Ll <0{l)Si, ||«.(t)||oo <0(l)<5i, z=l, 2. 

Remark 4.1. The estimate of ||wi;(i)||Li can also be obtained directly from (|3.14|) : indeed, since 

{pix - XiPi)it, L) = 
and the total variation of ut, l is bounded by , from (|3.12|l one gets 

/ \vix - Aiui|(s, L)ds < Si, 
Jo 

and hence ||wi(t)||Li < 2Si. 

To obtain the estimate on V2 from H3.14|l one has to start supposing 

|e(s, 0)\ds < Si. (4.7) 



With the same computations as before one gets ||u2(t)||Li < 0{l)Si. As it will be clear from the 
next sections, the assumption (|4.7(l actually leads to the estimate 

ft 

|e(s, 0)\ds < 0{l)6l, 



and therefore it is a posteriori well justified. 
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4.1.2 Integrability with respect to time 



The following lemma, which can be proved by a simple integration by parts, introduces a useful 
estimate we will widely use in the following. 

Lemma 4.1. Let P(x) be a non negative function defined on R and let q be a solution of 

qt + i>^q)x - qxx = s{t, x). 

Then the following estimate holds: 

4/ \q{t, x)\P{x)dx < [ \s{t,x)\P{x)dx+ [ \q{t, x)\{XP' + P"){x)dx 
dt In ./n ./n 



P'W)\ 



x=0 



Psign((7) {q^ - Xq){t) 



x=0 



Before applying the previous lemma, we recall that the boundary data of the scaled problem 
belongs to BV{0, +oo) and that the norms of u[,g and are bounded by Si. From the 
decomposition ut = wivi + W2r2, we immediately have 

\\wi{x = Q)\\l^o,+cx>) < Si \\wi{x = L)\\li(^o,+oo) < Si i = 1, 2. 

Moreover, in Section rOl we found that m; i 1, 2 can be decomposed as follows: 

wi = pix - Aipi + vix - Xivi 

W2 = P2x - X2P2 + V2x - X2V2 + e{t, x), 



(4.8) 



where the error term e(t, x) satisfies the estimate As we anticipated in Remark 14. II we will 

suppose 

|e(s, x)\ds <Si Vx e [0, L\. 

Since we will obtain an estimate of order S^ < Si, this assumption is a posteriori well justified. 
From the boundary condition H3.18|l {p2x — X2P2){t, 0) = and from the decomposition H4.8|l we 

get 



\V2X - X2V2\{s, 0)ds < 2Si. 



Similarly, one obtains that 



/ \vix - Aiwi|(s, L)ds < Si. 
Jo 

An application of Lemma [4. II with P = 1 and g = f 2 leads by observation H3.15|l to 



\v2x{s, L)\ds < 



t rL 



"'0 



\s2{s, x)\dxds + / \v2x ~ X2V2\{s, 0)ds + / \v2{0, x)\dx 



and similarly 



< Oil)Si + 25i + Oil)Si < 0{l)Si 



\vixis, 0)\ds<O{l)Si 



Let 2c be the separation speed defined by H1.5|l : the application of Lemma l4.1l with q{t, x) = V2(t, x) 
and 

1/c X < y 

P{x) = Py{x) = { y e [0, L[ 

exp I c[y ~ x)\ /c x > y 
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leads to the estimate 

\v2{s, y)\ds < i \v2{Q, x)\dx + - / / \si{s, x)\dsdx 
Jo c Jo Jo 

+ Py{Q) f \v2x - A2«2|(s, 0)ds + Py{L) [ \v2,is, L)\ds 
Jo Jo 

< + 0{l)Si < 0(l)<5i Vy e [0, L[. 

Analogously, we get 

t 

\vi{s, y)\ds < 0{l)6i yy e]Q, L]. 

The following proposition summarizes what we have proved so far: 

Proposition 4.2. Let Vi, pi i = 1, 2 be the solutions to the equations (ll.21() with the boundary 
conditions described in Section \H.'A Then it holds 

t pt 
\v2x- >^2V2\{s,Q)ds <25i, I \vi^, ~ \ivi\{s, L)ds < 5i, 

Jo 

\vi^{s,Q)\ds<0{l)8i, / \v2x{s,L)\ds<0{l)5i, 
Jo 

and 







\v,{s,y)\ds<0{l)5i, yy G [0, L] z = 1, 2. 
Further computations f Appendix IA.3. l)l ensure that 

\piAt, x)\ < 0(l)<5iexp(-ca;/2), |p2x(i, x)\ < exp (c(x - L)/2) . (4.9) 

The following proposition deals with other estimates of integrals with respect to time: the proof is 
quite long and requires the introduction of new convolution kernels. It can be found in the Appendix 

Proposition 4.3. In the same hypothesis of Proposition it holds 

I \v^x{s, y)\ds < 0{l)Si yy G [0, L] i = 1, 2 
Jo 

and ^ 

[ \w,is, y)\ds < 0{l)Si yye[0,L] ^ = 1, 2. 
Jo 

We also have 

ft 

\w^x{s,y)\ds<0{l)6l yye[0,L] i ^ I, 2. 



In the previous proposition the functions Wi are of course defined by relation ut — wifi + W2r2- 
Putting together Proposition 14.21 and 14.31 and the decomposition H4.8|l one gets 



/ \pi.-XiPi\{s,y)ds<0{l)Si, f \p2.~X2P2\{s,y)ds<0{l)5u Vy S [0, L], 
Jo Jo 



10 

and 



/ \pix - >^iPi\{s, 0)ds < mSi, / \p2x ~ >^2P2\{s, L)ds < mSi, 
Jo Jo 

where the constant m satisfies the hypothesis stated in Section 0] 
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The estimates obtained so far will be widely used in next sections and moreover allow to prove 
a bound of order 0{l)6f on some of the terms that appear on the right hand side of H3.13II : 



t pL 

Jo 



^{\Pi\ + \Pix\){\vj\ + \v.jx\ + \wj\ + \wjx\){s, x)dsdx 



(4.10) 



< 0{l)5i / (e-^^ + e^(-^)) / (1^,^-1 + \v,,\ + \w,\ + |u;,,|)(s, x)dsdx < Oil)S'f 



and 



/ / \Pix - AiPiKIp^I + \pix\){s, x)dxds 

JQ Jo 



< 0{l)di / e""^ + e"(^--^) / \pix - XiPi\{s, x)dsdx < 0{l)dl 



(4.11) 



4.2 Interaction functionals 

In this section we introduce three nonlinear functionals and we use them to bound those terms in the 
right hand side of due to interaction between waves of different families and those due to the 
fact that the speed cti is not constant. The form of the functionals is exactly the same considered 
in \7\ , with some more technicalities due to the presence of the boundary. 

4.2.1 Interaction among waves of different families 

We claim that the condition 



t i-L 



|Si(s, x)\dsdx < 0{l)Si 



t i-L 



\s2{s, x)\dsdx < 0(l)(5l 



implies 



I I + l^^^l) + l^^^^^lj ^)dsdx < 0{l)Sl 



We will prove only that 



Jo 



\viV2\{s, x)dsdx < 0{l)dl, 



(4.12) 



(4.13) 



because the other terms in H4.12|l can be dealt with analogously: see for example (H). 
Let 2c be the separation speed introduced in H1.5(l and let P{^) be defined as follows: 



l/2c 



One gets 
ds 



P{x - y)\vi{s, x)\ |w2(s, y)\dxdy ] < / \v2{s, y)\ 



Jo 

L 



P{x - y)s\gTLVi{vi^ - \ivi){s, x) 



dy 



+ 



|ui(s, x)\ 







P(x - y)signu2(w22; - A2U2)(s, y) 



dx 



y=o 



|w2(s, 



P'{x-y)\vi{s, x)\ 



x=0 
=L 



dy 



|ui(s, a:;)! 

JO 



y=L pL 

+ \viis, x)\ I P{x-y)\si{s,y)\dy 

y=0 JO 



P- 



P'{x-y)\v2{s, y)\ 
'{x - y)(^Xi{s, x) - A2(s, y)j + 2P"{x - y)j\vi{s, x)\ \v2{s, y)\dxdy. 
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One has 

P'(Ai - A2) + 2P" < 2(-cP' + P") = -5s=o, 0<P(s)<:^, 0<P'(s)<i 

Zc 2 

and moreover from the estimates of Proposition l4.ll and l4.2l it fohows that 

\vi^ - Xivi\{s, L) \v2{s, y)\dyds < 0{l)6f / - Aa^sKs, 0) / ^s, x)\dxds < 0{l)Sl 
Jo Jo Jo 

y)\ / |wi(s, x)\dxdyds < 0{l)dl : 
Jo 

this completes the proof of the estimate H4.13|l . 

With some technical computations, in Appendix I A . 3 . 3l it is proved 

/ / II + I^J-l) + < 0{l)Sl (4.14) 

which completes the proof of the estimate 

/ / H + + \w,\ + \w,,\^ + \w,\{\w,\ + x)dsdx < 0{l)Sl 

4.2.2 Length and area functionals 

To prove the estimate 

1^12:^1 — viWix\{s, x)dxds < 0(1)^^, 



"'0 



we introduce the curve 



and the related area functional 



Iff. , . . 1 '■^ 



-4(7)(s) = o / / \lx f\ lv\dxdy ^ - I I \vi{s, x)wi{s, y) - vi{s, y)wi{s, x)\dxdy. 

^ J Jy<x ^ Jo Jo 

The curve satisfies 

Ixt + (Al7a;)a: = ^xxx 

and moreover one has 

dAi^s) If f gjgj^|^^^(-g^ x)wi{s, y) - vi{s, y)wi{s, x)j ^wi(s, x)wi{s, y) - vi{s, y)wi{s, 

n(vi{s, x)wi{s, y) - vi{s, y)wi{s, x)j ^Ai(s, x)(vi{s, x)wi{s, y) - wi(s, y)wi{s, x) 
n(vi{s, x)wi{s, y) - ui(s, y)wi{s, x)^ [vi{s, x)wi{s, y) - vi{s, y)wi{s, x) 
sign(^wi(s, x)wi{s, y) - ui(s, y)wi{s, x)^ f Ai(s, y){vi{s, x)wi{s, y) - vi{s, y)wi{s, x) 



ds 2 Jq 

^ 2 



sig: 



'0 Jy 

-L fx 



sig: 



'0 "'0 

»L l-x 



1 

2 

1 

2 Jo Jo 
and hence 

dA \ , .... . . ^ . , 1 '■■^ 



-7- \ \vxy{s,L)wx{s,y)-v-Y[s,y)wYy{s,L)\dy- - \ \vxy[s, y)wx{s, y) - wiy(s, y)vi(s., y)\dy 



ds -2 .In ' ' ' ' ' ' ' ' " ^ 2 

r-L 

Ai(s, L)|t;i(s, L)ti;i(s, y) - ui(s, ?/)wi(s, 



1 

I , 1 rL 



vi{s, x)wix{s, x) - wiis, x)vixis, x)\dx + 2 / ri('5' ^)^ix{s, 0) - vix{s, 0)wi{s, x)\dx 
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Since ^(7)(0) < 0{l)6f, one obtains, using the estimates in Propositions 14.21 and 14.31 



t pL 



Jo 



dA 

\vi(s, x)wix(s, x) — vix{s, x)wi(s, x)\dx < — I ——ds + 0(l)Si < 0(1)5?. 
' ' ' ds 



The length functional of the curve (|4.15|) is defined as 



^{'l){s)= I \lx\dx^ I ^vl+wfdx, 



and wiU be used to prove the estimate 



"'0 



Wi 
Vl 



n 2 



Xdxds < 0{l)Sf, 



where x is the characteristic function of the set 

Wi 



(x) - XI 



< 3(5 



(see Section l3.1.3l for the definition of 6). 

We preliminary observe that the foUowing equahties hold: 



^^1 



^L'^r +'"1x^1 - '^VlxWixViWi ^ ^\jxx\ \lx\ - ilx, Ixx) 



\Mlx\x 



\lx 



{Xl^x, (Al7x)x) {ix, {Mlx)x) 



|Ai7x| 



17=. I 



{lx-,lxx)\ _{lx,lxx)'^ , {lx,lxxx) , \lxx\'^ 



From 7a;t + (Ai7a;)a; = ^xxx, one gets integrating by parts 



dL 

ds 



ilxxx, 7x) f {{Mix) a; J 7a;/ 



17a; 



17a; 



Hence, 
^ Jo Jo 



\Vl\ 



\lx \xx 



xdxds < 



< - 



{ix, Ixx 



\lxx? 



|7xP Jo \lx\ 

* |7a:a;P|7a:P " (7a;, Ixx)'^ 



\Mlx 



dxds 



\lx 



(4.16) 



dC 

ds 



ds + 







x—L 






f 


hx\x{s, X) 




f 


\Mlx\{s, X) 


Jo 




x=0 


Jo 





X—L 



ds < 0(l)(5i 



In the previous estimate we have used the fact that vi, wi, vix, wix are integrable with respect 
to time and that their integrals are bounded by 0{l)Si fPropositions 14.21 and l4.3|l . Since ||wi|loo is 
bounded by 0{l)Si, the previous estimate complete the proof of 14.16|l . 



4.3 Estimate on the error in choosing the speed 

The final estimate is the source term due to the cutoff function 9. Also this computation is similar 
to the one performed in 0, taking into account the fact that here we have a double boundary. In 
Appendix I A . 3 . 41 one can find the proof of the estimates 




(4.17) 
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This ends the proof of the estimate 

rt pL 



"'0 



\§,{s, x)\dsdx < 0{l)5l i = l,2, 



and hence of Theorem l4.1l 



5 Stability estimates 

In this section we prove the second part of Theorem ll.il completing the proof. Since the ideas are 
essentially the same as in the boundary free case, we will only sketch the line of the proof, paying 
more attention to the choice of the boundary conditions (which is the new element in this paper). 
The result of this section is thus: 

Theorem 5.1. There exist constants Li, L2 s.t. the following holds: let ^ be two solutions of 
the parabolic system 

ut + A{u)ux - ^ 0, (5.1) 
with initial and boundary data mJ, u^q, and Uq, u^q, u^^ respectively. Then 

\\u^(t) - -u^(s)||ii(o,L) < Li{\\ul - Uo||li(o,l) + - Ubol|Li(o, +00) + \\uIq - uLlUi(o, +00)) 
■L2{\t~s\ + \Vt-^s\). 

(5.2) 



5.1 Stability with respect to initial and boundary data 

We will prove that, in the hypothesis of Theorem 15. II 

\\u^{t) - U^{t)\\L^n.,L) < Li(\\ul - wg||Li(o,L) + 11^60 - uIo\\l^O.+oo) + \\ulo - U6lI1l1(0,+oo)) (5.3) 

Let z{t, x) be a first order perturbation of a solution u{t, x) of (|5.1|l . By straightforward computa- 
tions one gets that z satisfies 



zt + {A{u)z)^ - Zxx = {DA{u)ux)z - {DA{u)z)^ 



(5.4) 



To prove Theorem 15.11 it is enough to prove that any first order perturbation z{t, x) satisfies the 
bound 

lk(i)llLi(o,i) ^ Li(^\\z{t = 0)||li(o,l) + \\z{x = O)||ii(o,+oo) + \\z{x = i)||Li(o, +00)) • (5-5) 

Indeed, provided l|5.5|l holds, a homotopy argument which can be found in [SJ El gives then the 
Lipschitz estimate (|5.3ll . 

To prove H5.5|l it is convenient to introduce the auxiliary variable 



which satisfies the equation 



Tt + {A{u)T)x - Txx = 



T = - A{u)z, 



DA{u){ux^z — z^Ux) — A{u) 



DA{u)(^Ux ® z — z ® Ux) 



+ DA{u) {ux (g) T) - DA{u) (ut (g) z) . 



(5.6) 



Let 20(2;), zi,o{t) and Zf, ^{t) be the initial and boundary conditions we impose on z: since the final 
goal is to apply (|5.5|) in the homotopy argument, it is not restrictive to suppose that zq{x), zi,o{t) 
and zi,L{t) satisfy the same regularity hypothesis as u. Indeed, the solution z of l|5.4|l that is used 



29 



in the homotopy argument is on the boundaries and at t = just the difference of the solutions 
and of H5.1|l . 

Hence we will suppose that zo{x), Zbo{t) and ztLit) are regular and that d'^zo/dx'^, d'^z^o/dt'^ 
and d'^Zhh/dt'^, k = 1, . . . n are integrable and have a small norm. Moreover, if — UQ||ii(o,L), 
\\u\f^ — 1*^011^1(0, +cx;) or \\u\j^ — u^^||ii(o, +oo) a^c infinite, then (|5.3(l holds trivially, and therefore 
we can suppose that zq € iv^(0, L), ZbO: ZbL & ^^(0, +oo). 

From the hypothesis on zq it immediately follows that T(i = 0) is regular and small in and 
sup norm. 

As in the proof of the BV bounds on the solution u, the crucial step to show 1)5.5(1 is the 
introduction of a suitable decomposition along travelling waves and double boundary layers: note, 
moreover, that Ux satisfies equation ((5.4(1 . Hence, it seems promising to decompose z along the 
same vectors ri{u, Vi, Oi) and fi{u, pt) that appear in the decomposition 1(3. lUI) of u^- This choice 
actually leads to non integrable source terms. We will therefore allow the vectors employed in the 
decomposition of z to depend not only on the solution u, but also on the perturbation z itself: 

z = zifi{u, vi, Ti) + Z2r2 + qifi{u, pi) + g2?'2 

T = Wl, Tl) + L2r2. 

In the previous expression the speed of the travelling waves described by the vector fi is not ai , but 
The function 6 is the cutoff 



s 



if\s\<S 

9{s) = <( if |s| > 3(5 S< 



smooth connection if S < s < 3S 
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The proof of 1(5.5(1 is from now on very similar to that of the BV bounds: one inserts the previous 
decomposition in the equations 1(5.4(1 and ((5.6(1 and obtains the equations: 

Zu + (AiZi)j, - Zixx = Z2t + {X2Z2)x - Z2xx = Si(t, x) 

qit + iXiqi)x - qixx ^ q2t + C^2q2)x - q2xx = (5-7) 

Lit + (Alil)i, - Llxx = S 3(i, x) i2t + (A2t2)a; " l^2xx = Sjl^, x) 

As in the proof of the BV bounds, to prove ((5.5() it is sufficient to show that the condition 

lk(5)||Li(o,L) <C<5i Vse[0, 

implies 

/ / \si{s,x)\dxds<0{l)5l i = l, 2, 3 
Ja Jq 

and suitable bounds on the boundary terms. Moreover, in the proof of the previous implication it 
is not restrictive to assume 



[ [ \Si{s,x)\dxds<0{l)Si i = l, 2, 3, 
Jo Jo 



because a posteriori one finds a bound of order 6f. 

Actually, one could observe that while the equations for Ux and Ut have no source term (see 
Appendix I A . 2 . II for details), the equations ((5.411 and 115.6(1 have nontrivial source terms. However, 
one can show that both the source terms in ((5.4() and 1(5.6(1 and the other terms that contribute to 
Sj, i = 1, 2, 3 can be bounded by an expression analogous to the one that appears on the right side 
of 1(3. 13() . The computations that ensure such an estimate are quite similar to those performed in 
the proof of Section ITTI 

The proof of ((5.5() can therefore be completed with the same tools described in Paragraph 01 
hence we will skip all the details. 
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5.2 Stability with respect to time 

Let u{t, x) be a solution of (|5.1|) : from Proposition 12 . 21 and the observations that foUow one gets 



Let ti < t2- the estimate above imphes 



\\u{ti) - M(t2)||Li(0,L) < 



du 
'dt 

■t; 



0(1)^1 t>l. 



dt< [ {0{l)\\u^(t,x)\\Li + \\u.,x{t,x))\\Li)dt 



<0{l) f\6i+Si/Vi)dt<0{l)Si\ti-t2\+0{l)Si\Vh-Vh\ ^^-^^ 
Jtl 

<L2{\h^t2\ + \Vh-Vh\)- 
This completes the proof of Theorem 15 . II and hence of Theorem ll.il 

6 The vanishing viscosity Umit 

In this section we prove Theorem ll.2l The proof proceeds in two steps: first, by using the results of 
Theorem ll.il we obtain that there exists a subsequence of solutions u'^ to the problem 

ut + A{u)ua; ^ 0, X e]0, l[ t e]0,+oo[ 

u(0, x) — uq{x) 

u(i,0) = -Obo(i) u{t,l) = Ubi{t) 

which converges to a Lipschitz semigroup. Then we use the machinery of viscosity solutions to 
complete the proof, showing the uniqueness of the limit. In particular, we exhibit explicitly the 
boundary Riemann solver. 

Let [mo, UbOi Ubi] the solution of the system H1.3|l : from Theorem 14.11 one gets that the total 
variation of the solution of system (|1.16|) is uniformly bounded with respect to time and hence, by 
a change of variables, [uq, UbO: ^bi] satisfies 

Tot Var{p/ [mo, "6 0, "6i]}, [uo, Ubo, um] (a;) < 0{\)5i > 0, a; e [0, Z], e > 

and for any uq G Uq Ubo, Ubi <EUb- By Holly's theorem, for every sequence £„ 0+ and for any 
i > there exists a subsequence, which we still call £„ for simplicity, such that p/"['Uo, UbOj Ubi\ 
converges in L^(0, I). The stability with respect to time and to initial and boundary data ensures 
that, by a standard diagonalization procedure, one can find a function 

p : [0, +oo[xUi^ xUb y-Ub Vq 

{t, Mo, Ubo, Ubi) i-> Pt[uo, Ubo, Ubi] 

such that, up to subsequences, 

p/"(i)[uo, Ubo, Ubi] Pt[uQ, Ubo, Ubi] L^(0, I) yt > 0, uo e Uq, Ubo, Ubi G Ub- 
Moreover, one can verify that the function 

S : [0, +oo[ X Z^o X X VaxUbxHb 



{t, Uo, UbO, Ubl) 



Pt[Mo, UbO, Ubl], Uo{- +t), Ubo{- +t), Ubl{- +t) 



(6.1) 
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satisfies the semigroup properties, together with the Lipschitz estimate 



Pt[uo, Ubo, Ubi\ -Ps[vo, Vbo, Vbi\ 



<Li \\vo - UoIIlI(OJ) + \\vbO - Wf,o||Li(0,+oo) 



(6.2) 



\Vbl - Uw||l1(0,+oo) j + L2\t - s\, 

We now make use of the tool of viscosity sohition, which was first introduced in 

6.1 The Riemann solver and the boundary Riemann solver 

A crucial step in the proof of the uniqueness of the vanishing viscosity limit is the local description 
of the vanishing viscosity solution in case of piecewise constant data, which however has an interest 
in its own. The aim of this section is to characterize the limit as e„ 0+ of the solution of 



Ut + A{u)Ux = SnUxx 

u+ X > 
a; < 
u{t, I) 



u{Q, x) 
u{t, 0) 



(6.3) 



"6 



Ubl 



where , u~ , Ubo and Ubi are constants. In the following, we will write "solution to the Riemann 
problem" meaning "vanishing viscosity solution to the Riemann problem". 

In El it is shown that the solution of 1)6.3(1 is defined locally: to solve 1)6. 3|l it is therefore 
sufficient to characterize the vanishing viscous solutions in the following three cases: 



1. the Cauchy problem with datum 



uq{x) 



X <0 
X > 



2. the boundary problem at x = 

r u(0, x) = uq 

\ U{t, 0) EE UbO 

3. the boundary problem aX x — I 

r U(0, X) EE Uq 
\ u{t, l) = Ubl 

The second and the third case are clearly analogous, and therefore in Section IfT^ we will deal only 
with the second one. In the following section, instead, we will recall for completeness the essential 
steps of the construction of the solution in case 1: we refer to T for an exhaustive account. 

In any of the three cases the crucial step is the definition of two families of admissible states, as 
it will be clearer in the following. 



6.2 The non conservative Riemann solver 

Since in this case the construction of the first and the second curve of admissible states is the same, 
we will describe only the construction of the first curve T^u^ of the states that can be connected by 
waves of the first family to a right state Ur- For a general reference, see [7]. 
Consider the family T C C°( [0, s]; R" x R x R) of curves 

T ^ {u{t), Vi{t), Cri(T)), T e [0, S] 

with 

\u(T)-u*\<e, |«i|<e, |At-ai(T)| <£. 
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The function /i (r) related to the curve 7 G T is defined as 



Let fi be the generahzed eigenvector of the travelling waves of first family (see Section r?.1.2l for the 
proper definition of fi ) . By the contraction map principle, one can show that if s is small enough 
then for any r € [0, s] there is a solution (u, vi, cti) of the following system: 



u{t) =Ur+ ri{ui<;), vi{<^), ^ii<;)d<; 





vi(t) = conC[o,s]/i(T) - /i(t) 
^ (iconc[o,s]/i 
^^(") = 

We indicate with conc[o,s]/i the concave envelope of /i in the interval [0, s]. 

The curve of admissible states passing through Ur is defined as T^Ur = u{s). Indeed, let 

T^Ur x/t < ai{s) 
u{x/t) = ^ u{t) ai{T) = x/t 
x/t > ai{0) : 

one can show that any sequence of vanishing viscosity solution of the Riemann problem with data 
{ur, T^Ur) converges to u. Moreover, the curve T^Ur is Lipchitz continuous. 

6.3 The boundary Riemann solver 

In this paragraph we will construct the vanishing viscosity solution in case 2. We will proceed as 
follows: we will construct two curves of admissible states and and given a right state uq and 
a left state w^Oj we will show that there is a couple (si, S2) such that 

^ii ° zI^uq = ub(). 

The waves of the second family are entering the domain: it is therefore quite reasonable to suppose 
that they are not influenced by the presence of the boundary and therefore the second admissible 
curve will be the one defined in the previous paragraph, Z'^uq — T^uq. Let u — Z^^uq be the value 
reached throughout the waves of the first family. 

The waves of the first family are leaving the domain and are therefore affected by the boundary 
datum. To understand their behavior, it is convenient to focus the attention on the boundary layers 
of the first family, i.e. on the solution of 

Uxx = A{u)Ua; (6.4) 

that are exponentially decreasing to an equilibrium as x — s- +00. One can now go back to the 
problem 

and let e 0+. Since u^{x) = u{x/e), we get 

lim it^(O^) = lim u{x). (6-5) 

Such a behavior is illustrated in figure 01 

The value lim(,^Q+ u^(0+) is the state reached throughout the waves of the second family: we 
called it u. It also represents the trace of the hyperbolic limit on the boundary x = 0. From ]<j.b\i 
it follows that the states which can be connected to u by boundary layers are the initial points of 
orbits that decrease exponentially to u, i.e. that lay on the stable manifold throughout u. 
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The stable manifold at the equilibrium point (u, 0) of the system 

Px = A{u)p 

is parameterized by the projection pi of p on the stable space. Passages analogous to those in Section 
I3.1l ensure that the stable manifold is characterized by the relation p — pifi{pi) for a suitable vector 
function 

1 

n 



/(Pi) 

One imposes wi(+oo) = (/i, u) and from the second equation gets 



wi(0) = {h, u) -pi(0)exp ^ j \i{ui{pi{{)), x)^dx 



Since Ai < — c < 0, the previous map is invertible and one can express pi(0) as a function of ui(0). 
The inverse map is clearly regular. 

We parameterize the stable manifold by si := ui— < h, u > and obtain (for some suitable 
regular function z) the map 

^Mt^O> (6.7) 



z(si) 

defined on a small enough interval [0, s]. 
The vanishing viscosity solution of 

ut + A{u)ux = 

(6.8) 

u{t, 0) = Uq u(0, x) = UbQ 

can be constructed patching together the curve described so far. Let 

Mfeo = o Tl_uo ■■ 

thanks to a version of the implicit function theorem valid for Lipschitz maps (see |18|), one can 
reconstruct from uq and uto the couple (si, 52)- The vanishing viscosity solution of H6.8() is then 
given by 

T^^uo x/t < 0-2(52) 
u{t, x) — u{t) cr2(T) = x/t 

Uq x/t > 0-2(0). 
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One gets in particular that the trace of the solution at x = is not necessarily the boundary value 
Ub Oi but it is the intermediate state T^^uq. 

Remark 6.1. In the case of systems in conservation form, with only linearly degenerate or genuinely 
non linear fields, a boundary Riemman solver was introduced in |21j . In that paper, it was introduced 
the following admissibility condition on the trace u{t, 0+) = u oi the solution of 1)6.8(1 : the solution 
in the sense of Lax |26) of the Riemann problem 

Ut + I{u)x = 

Ub a; < 



m(0, x) 



u a; > 



is composed only of waves with non positive speed. Such a condition is in general different from 
(|6.7(l and therefore the two boundary Riemann solvers do not coincide. 

On the other side, in 25, 31, it was considered a quite general boundary condition: more 
precisely, let N be the dimension of the system 

ut + f{u)x = 

and let p be the number of positive eigenvalues of Df(u), which is supposed to be constant. Let 
b : M.-^ ^ RP be a regular enough function such that Dh{u) is injective on the space generated by 
the p eigenvectors of Df{u) associated to positive eigenvalues; then, given g : [0, +oo[^ W, the 
boundary condition considered in [231210121 is g{t) = b[u{t, 0+)). Such a definition, which in the 
original papers was introduced in the case of conservative systems with only linearly degenerate or 
genuinely non linear fields, is compatible with the boundary Riemann solver defined by the vanishing 
viscosity limit. Indeed, in our case N = 2, p = 1: let b{u) be equal to the coordinate of u along 
the curve of admissible states T^uq, i.e. let b{u) = S2 H u = Zl_^ o T^^uq. Moreover, let g be the 
coordinate of Ub along the same curve: with this choice, the condition 

u{t, 0+) = T^^uq Ub = Z\ o T^^uo 

is equivalent to g{t) — b(u{t, 0+)). 

6.4 Viscosity solutions 

Before giving the definition of viscosity solution we have to introduce some preliminary notation; 
moreover, in the following we will use the spaces Uq, Ub, 2?o that have been defined in the introduction 
(equation Hl.l()|l and previous lines). 

Let u{t, x) be a function such that, for any t, u{t) € T>o: given a point (r, e]0, ^[x[0, +oo[, 
let A'^ = A(u{t, 1^)) and let U^^ ^ be the solution of the linear Cauchy problem 

Wt + A'^Wx — w(0, a;) ~ uir, x). 
Viceversa, let J/'^ ^ be the solution (defined in Section IH?^ of the Riemann problem 

Ut + A{u)ux = 

in \ f "('^: C") a; < 
m(0, a;) = ; ^ , ( ^ 
^ ' ' \ u{t, ^+) a; > 

The previous limits are well defined, since u{t) G BV{0, I). Given a function Ubo G Ub, the 
definition of uf can be extended naturally to the case £ — 0: it is enough to define uf - <, 
as the solution (described in Section l6. 3(1 of the boundary Riemann problem 

Ut + A{u)ux = 

m(0, a;) EE u(r, 0+) u{t, 0) = Ubo{T+). 
Given a function Ubi G Ub, the definition of U^^ -^^ is clearly analogous. 
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Definition 6.1. Let u{t, x) such that for any u{t) £ T>q and such that the function 1 1-^ u{t, ■ ) is 
continuous in Lj^^ and let ui,q, Ubi € Ub and uq £Uq. 
Then u is a viscosity solution of the system 

Ut+A{u)ux^Q, a; e]0, < e]0, +oo[ 

u{fd, x) ^uo{x) (6.9) 

u{t,0) ^ Ubo{t) u{t,l) = Ubl{t) 

if and only if the foUowings hold: 

(i) u(0) = Mo 

(ii) for every /3 > and for every point (r, ^) with ^ 7^ 0, Z 

lim - / \u{t + h, x) - Tj} ^ ^Ah, x - S,)\dx ^ Q 

(iii) for every /3 > and for every r > 



1 



min{/, 



lim — / \u(t + h, x) — U, - -.{h, x)\dx = Q\ 



and 

lim — / \u(t + h, x) — Uf - -.(h, x)\dx — 

/i^0+ ft Jniax{0J-/3?i} ^ 

(iv) there exist constants C and (3' such that for every point (r, ^) with ^ 7^ 0, Z and for every p > 
small enough 



1 



min{/,^+p-,3'/i} 



limsup^ / \u{T + h,x)-U^^^^^{h,x-^)\dx<C[TotYsiT{u{T),]^-p,C + p[)] . 

The previous definition may appear a bit complex: note, however, that, since p and h can be 
arbitrarily small, it is not restrictive to suppose 

max{0, ^-ph} = ^-f3h min{Z, ^ + f3h} = ^ + f3h 

maxjo, ^-p + P'h} ^^^p + f3'h min{l, ^ + p- (3'h} = ^ + p-f3'h 

maxjo, I - h(3} = I - h(3 min{?, hp} = h(3 

The definition of viscosity solution ensures, roughly speaking, that a function is well approximated 
by the solution of a suitable linear problem and of a suitable Riemann problem. 

The following proposition ensures that viscosity solutions coincide indeed with vanishing viscosity 
limits. The proof is very simile to that of the analogous property stated in 7 (Lemma 15.2, page 
308) and will be therefore omitted. 

Proposition 6.1. Let uq € Uq and Ubo, Ubi G Lib- Let pt{uo, Ubo, Ubi) be a vanishing viscosity 
solution of the system (|6.9() ; then pt{uo, Ubo, Ubi) is a viscosity solution of the same system. 
Viceversa, if u(t, x) is a viscosity solution of the problem (|6.9() then 

u{t) = pt{uo, Ubo, Ubi) \ft>0. 

From the previous result it immediately follows the uniqueness of the semigroup: indeed, let 
by contradiction P((?2o, Ubo, Ubi) andp(('Uo, Ubo, Ubi) be two different vanishing viscosity solutions. 
The function pI{uq, Ubo, Ubi) is hence a viscosity solution of problem by the first part of 

Proposition 16. II Then p]{ua, Ubo, Ubi) = p1{uq, Ubo, Ubi) for any i > by the second part of the 
proposition. 
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A Appendix 

A.l Appendix to Section [21 
A. 1.1 Proof of Proposition [!i.ll 

In the following, for simplicity we will suppose A* = A2 > 0, since the case A* = A* < is analogous. 
We denote by 

r^^* it, X, v) = (1 - e-'^y")G{t, x-y- X^t) 

the solution of the equation 



Zt + A2ZJ: — Zxx — 



(A.l) 



in the first quadrant with zero boundary datum and Cauchy datum 6y. The following estimates 
have been proved in [7]: 



Wr^Ht, y)llLHo,+oo) <o(i) 
l|r^'(i, y)\\LHo.+oo)<^ 



4-oc 



r>{t, X, 



dx 
dx 



< 0(1) VteM+ 



(A.2) 



Let G{t, x) = exp(— a; /4t)/2v7rt: we will use the notation 



G{t, x~X*2t) = G^^{t, x). 

The estimate on the norm of A^'^ in Proposition l2.1l can be obtained via the maximum principle 
applied to equation l|A.l(l : indeed, 



0<A^'^{t,x,y)<G-^^t,x-y) Vi > y e]0, L[, 



and therefore \\A^^{t, y)\\Li < 1 



To prove the estimate on the norm of A^^ it is convenient to write A'''^ as 

A^2 (t^ X, y) =r^2 (t,x,y) + (t) ^'^ {t, x, y) ^ ^G{t, x + 2mL -y)-G{t,x + 2mL + y) 



with 



4>^'^{t, X, y) = exp I ^(x-y) 



AS 



(A^) 



*\2 



Since A2 > 0, for m > it holds 
^ (^0^^ (t, X, y)G{t, x~y + 2mL)^ 
= (l)^Ht, X, y) 



-^G'(i, x-y + 2mL) + Gx{t, x-y + 2mL) 



< \G7 {t, X - y + 2mL)\, 



and similarly 
dx 



< \GxHt: x + y + 2mL)\. 



exp (^^(a; - y) ~ G{t, x + y + 2mL)^ 

The terms of the series corresponding to to < can be estimated as follows: let n -.^ —to then 
X, y)G{t, x-y- 2nL) 



d_ 

dx 



^cj,^Ht, X, y) 



Gx{t, 2nL-x + y)- ^G{t, 2nL - x + y) + \^G{t, 2nL - x + y) 



< 



{t, 2nL-x + y) + A^ G^^ (^t, 2nL-x + y) 
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and similarly 



d 

— ( (j){t, X - y)Gx:{t, x + y- 2nL) 



Since ||G^' < 0(l)/%/t, one obtains 



< 



G^^ (i, 2nL-x-y) + G^^ {t, 2nL ^ x - y) 



< 



(t, X, y)\dx + 



2L 



G^Ht, z + y)\ + \Gi^t, z~y)\ ]dz 



A* 



Gt'{t, z + y)\ + \Gi'^t,z-y)\ ]dz + X 



\GPit, z + y)\ + \Gi^t,z-y)\ ]dz < 



Vt 



In the following estimates, we will suppose y < L/2: by symmetry this is not restrictive. Observe 
that, for y < L/2 

x + y-2L<-L/2<0 Va; e [0, L]. (A.3) 

This assumption corresponds to the fact that the most singular part in A"^ is collected in F'*' , i.e. 
it is given by G{t, x + y) — G{t, x — y). If y > L/2, then the most singular part would be given by 

G{x~y)~G{x + y-2L). 
One has 

it^ x,y)= f r^* (i, X, + [ Mt, X, C) Yl ^ - ^ + 2"^i) - G(i, x + £, + 2mL)\ d^ 
Jy Jy Ln 



+ I (f>{t, X, ^) \G^{t, x-^ + 2mL) ~ G:,{t, X + + 2mL) d^ 

+ f U^it, X, c) ^ - ^ + 2™i) + Ht, X, Y 2; - e + 2mi)|de 

: / {t, X, £.)d^ + Y <^(^' ^' 2/)G'(i, a; + y + 2mL) - Y ^i^^ L)G{t, x + L + 2mL) 

•^y m^O m/0 

- I Y ^2^^^^ ^' OG{t, x + i + 2mL)di + Y ^(''^^ ^' y)Git, x-y + 2mL) 

•'y m^O m.^n 



- Y '?^(*' ^' L)G{t, X-L + 2mL). 

The integrability of the first term follows from (|A.2|I . the other terms are clearly integrable because 
of the quadratic exponential decay of the heat kernel G: hence || A'*'^ (i, y)||ii < 0(1). 
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The function A^^ can be written as follows: 



A^' {t, X, y) = / r^l{t, X, + -fHt, ^, y)G{t, ^ + y + '^mL) 



AS 



^ (t){t, X, y)Gx{t, x + y + 2mL) - ^ y'^^C^' 2;, L)G{t, x + L + 2mL) 

(j>{t, X, L)G^{t, a; + L + 2mL) - Y 2~^^^' ^' ^ + ^ + 2mL)d^ 

L 

X, —y + 2mL) 

Y ^{i^ X, y)Gx[t, X, -y + 2mL) - ^ -^(l>{t, x, L)G{t, x- L + 2mL) 



- Y '/'(^i X, L)Gx{t, x- L + 2mL), 
and hence with computations similar to those performed before one gets 

\\A^J{t,y)hi<^ yt<l ye]0,L[. 
If one derives the explicit formula (|2.6|) for J'*' ^ and then integrate by parts gets 



tA: l 



x)\dx = 



A*Ce^^^da; - A^C / A''^ {t, x, y)e^^ydy 



dx < 



where we have used the estimate ||A'^2||^i < 0{l). By symmetry it follows || < 0{1) 



Deriving Jx^ ^ one obtains 



Jxl'^WIUi< / \G{X;fe^'-^\dx + CX. 



* \ 2 ^XXx I 



A^Ht, X, y)e^'-ydy 



dx < 



0(1) 

Vi ' 



thanks to the estimate on ||Aa;^||. By symmetry one gets \\Jxx < 0{l)/Vi: this concludes 

the proof of Proposition 12. II 



A. 2 Appendix to Section |^ 
A. 2.1 Explicit source terms 

We want to find the equations satisfied by the quantities vi, V2, Pi, P2, wi, W2' we will use the 
decomposition 



Ux = vihiu, vi, (Ti) + V2r2 +piri{u,pi) +^2*^2 

Ut = Wifi{u, Ui, (Ji) + W2r2 

and insert it in the parabolic equation 

Ut + A{u)Ux - Uxx = 0. 

A derivation w.r.t. a; gives 

fix = Dri (Uifi + t>2r2 + Plfi + P2f2) + Vixfiy + CTixfio 
fix = Df 1 [Vifi + V2r2 + Plfl + P2f2) + Plxflp- 



Vi 
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Recalling that 

A2 := A2 -^1(^2, Drir2) 

A{u)fi = wiDfifi + Aifi + vi{Xi - ai)fiy 

A{u)fi = piDfifi + Aifi +piAifip 

one gets 

Ut = - A{u)ux 

= vixri + vifix + Pixfi + pifix + V2xr2 + P2xr2 - viA{u)fi - piA{u)fi - \2V2r2 - \2P2r2 
= {vix - Al^;l)(fl + vifiy) + vfaifiy + viV2Bfir2 + wipiDfifi + vip2Bfir2 + viaixfia- 

+ {Plx - AlPl)(fi + Plfip) + ViPiBfifi + V2Pl'Dfir2 + {V2x - >^2V2)r2 + {P2x - MP2)r2- 

(A.4) 

We multiply the previous expressions by ti and by i?2, the vectors of the dual basis of (fi, r2): we 
obtain 

Wi ^ Vix - XlVi + Plx - \iPi 

W2 = V2x - A2W2 + P2x - A2P2 + e(t, X), 

where the error term e{t, x) satisfies the estimate (|3.13|) in Paragraph 8.2.2. 
Deriving (|A.4p . one obtains 

Utx ^{vixx - (AlWl)a:)ri + {vi[vixx " (AlUl)a;) + 2vix{vix " AiUi) + {vlai)x^fiy 

+ {vi{vix ~ AiUi)^Dfifi + {v2{vix - AlWi) + {V1V2) ^T)fir2 

+ {pi{vix - AlWl) + {Vipi)x^^fifi + {p2{vix - AlWi) + {Vip2)x}j^fir2 

+ {crix{vix - AlWi) + {viaix)^fia + {vi[vix - AlWi) +vla-i^{fiy)x 

+ WlU2(Dfir2)a; + ViPi{T)fifi)x + ViP2{T)fir2)x + Vi(Tix{fla)x (A. 6) 

+ [pixx - (AlPl)a;^fi + {pi{pixx - [XlPljx) + '2pix{pix - AlPl)^fip 

+ {vi{pix - Alpi) + {vipi)x^T)fifi + {v2{pix - AlPi) + (w2Pi)x)Drir2 

+ {pi{pix - Aipi)^Dfifi + {p2{pix - Aipi)^Dfir2 + {pi{pix - Aipi)^ {fip)x 

+ ViPi{Y)fifi)x +W2Pl(Drir2) + {v2xx - {Mv2)x^r2 + {p2xx - (A2P2)x)?'2- 

On the other hand, 

Uxt = vitfi + vifu + V2tr2 +Pitfi +Pifit +P2tr2, (A. 7) 

where 



fit = Bfiiwih + W2r2) + vitriy + (Jufia 
fit = Dfi(wifi + W2r2) +Pitfip. 



(A.8) 
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We equal (|A.6|I and HA.7|) and we use HA.8|) . obtaining 

= Utx~ Uxt 

+ {vi(vi^ - Xivi) - uiWi^Dfifi + {v2{vix - Aiwi) + (^1^2)0; - viW2^T)fir2 

+ {pi{vix - Aiwi) + {vipij^Vififi + (^P2(wij: - Ai-wi) + {viP2)j^^fir2 

+ - Ait;i) + {viaix)^ - auv-^fia + {vi^vi + -^^(o-i - Ai)^ {fiv)x + wiD2(Dfir2)a: 

+ Vipi{V)fifi)x + WiP2(Dfir2)a; + wicri2:(ricr):c + [pixx - {XiPi)x -Pitjh 

+ (pi{Plxx - iXlPl)x -Pit) + '2pix{pix - AlPl)^fip + (vi{pia: - XlPl) + (wiPl)a; " WlPl^Dfifi 

+ (v2{pix - Aipi) + {v2Pi)x ~ W2Pij'Dfir2 +Pi{pix - Aipi^Dfifi + (^P2{pix - Aipi)^Dfir2 

+ (piiPlx - Xipi)j {fip)x + ViPi(Dfifi)x + V2Pl(Dfir2)x + (v2xx - iX2V2)x ~ V2t^r2 
+ (P2xx - {X2P2)x -P2tjr2 

= (yixx - (Aiwi)a; - vit^fi + (^pixx - iXiPi)x -pitjh 

+ (v2xx - iX2V2)x - V2t^r2 + {p2xx ~ {X2P2)x " P2t) f 2 + Si{t,x). 

We can therefore impose the conditions 

vit + {Xivi)x - vir^.^, = 

Pit + {XlPl)x -Plxx = 0^ 

V2t + iX2V2)x - V2xx = (^2 , 2;) , Si (i , x)) =Si{t,x) 
P2t + (X2P2)x - P2xx = 0, 

where (^1,^2) is the dual basis of (fi,r2). 

To compute the equations satisfied by wi, W2 we will use 

utt = u^xt - {A{u)ux)^ 

= utxx ~ {A{u)ut)^ + DA{u){uj, ®ut-ut® Ux), 

which follows from 

{A{u)ux)^ = DA{u){ut (E) Ux) + A{u)uxt 
[A{u)ut)^^ = DA{u){ux <E) ut) + A{u)utx- 

Straightforward computations ensures that 

Uxt - A{u)ut ^{wix - Aiwi)fi + wi{vix - Aiwi)fi„ 

+ WiViaifiy + wiV2'Dfir2 + WipiDfifi + wip2T)rir2 

+ WiCJixfla + iw2x ~ \2W2)r2 

and 

DA{u) [ux®ut — Ut®Ux] — viW2fi (8) r2 + V2Wir2 ® fi + PiWifi ® fi + piW2fi ® r2 + P2Wir2 ® n 



WiV2ri ®r2— wipiri (g) ri — WiP2ri ® r2 — W2Wi?'2 ® ri — W2Pir2 ® ri. 

(A.9) 
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Hence 

Utt = Witfl + W2tr2 

= -wlBfifi - wiW2Dfir2 - wivitriv - wiCFufic, + [wixx - [Mwi)3^ri + {vi{wix - AiWi)^Dfifi 

+ {v2{wi.j, ~ \iWi)^Bfir2 + {pi{wix - AiWi)^Dfifi + {p2{wix - AiWi)^Dfir2 + {vi^iwix - \iwi)^fi^ 
+ {^^ix{wix - AiWi)^fiCT + {wix{vix - Ait)i)^fi„ + {wi{vix - Aiwi)^ {fiv)x + wiviai{fiy)x 

+ {wiVi(7i)xfiv + {wiV2)x^fir2 + WiV2{Iifir2)x + {wipi)xT)fifi + WiPi{Dfifi)x 

+ {wiP2)x^fir2 + WiP2{^fir2)x + {Wiaix)xfla + Wiaix{fla)x + iw2xx " (A2W2)a:)''2 

+ {PlW2)x^fir2 +PlW2{T)fir2)x 

+ DA[u) {viW2fi (g) r2 + V2Wir2 ri + piWifi (g) fi + piW2fi ^ r2 + P2Wir2 fi 
- WiV2ri ®r2- wipifi (g) fi - wip2fi ® r2 ~ W2Vir2 <Si ri - W2Pif2 ® ri 

= - (AlWl)^:)^! + (w2xx - (A2W2)x)?'2 + S2(t, 2;). 

One can check that, since A is triangular, 

{£i, DA{u){ux (g) ut - ut (E) Ux)) =0 
and therefore the equations satisfied by Wi i ~ 1, 2 are 

Wit + (AiWi):e - Wlxx = 

W2t + {>^2W2)x - W2xx = {hit, x), S2(t, x)) = hit, x) . 

A.2.2 Proof of Proposition lOI 

Equation 1)3.10(1 and 1)3.11(1 ensure that, since, 



(A.IO) 



then 



Vl 



\vlaix\ = ©(l)!^!:^^! - VixWi\, 



Clx ^ 



Vl 



< 36. 



Most of the terms in Si(t, a;) i = 1, 2 and e(t, a;) can be directly reduced to those in ProDOsition l3.ll 
The terms which requires some technicalities are: 

1. 

\Plx - AlPlI |(^2(U, Vl, tTi), pi))\ < C(l)(|pi| + I^iDIpIo; - AiPi|. 

Indeed, 

lihiu, Vl, ai), hiu, pi))\ < m, h - rl)\ + \{i2 - 1*2, rl)\ < 0(l)(bi| + 

We have denoted by r* the first eigenvector of the matrix and by (^i, £3) the dual base 

of irl, r2). 
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2. 



< 



"ivix I wi - {Pix - Aipi) + 2vivia;ai + vfaix 



0{1) 



2uia:(AlPl -pix) 



VlxWi - ViWu 



3. 



(^ix{vix — Ai^i) + {viaix)x — critVi \ria\- some computations ensures that 

( — ) {vix - Mvi) + Vix(—] +Vi(—] -[ — ]vi=0. 
\Vi/x \Vi / x \Vi/xx \Vi / t 

Hence, since Iri^l = C(l)|ui|, one gets 



(Tlx{vix - AlUl) + {Viaia:)x - (TlfVl \fia\ < 0{l)Xnx, 



4|A5:-u,i/vi|<3<5}"l 



(-) 



4. I - Wiait + (Tixiwia; - XlWi) + (Wicrix)xl \fla\ ■ 

since 

-.v,9'(^] +w,xe'(^) -Mw,0'(^) +w,.9'(^) +n.,0'(^ 

\VlJt \VlJx \'"^Jx \^lJx V^^l 

one is left to the estimate 



Vi 



-w-l/v-l\<35} 



5. \viaix{fia)x\ '■ first of all, we observe that that 0'{s) ^ implies < and therefore 

WlxVl - VixWi 



IVK^ix] 



\v.0'\ (^) 



\vie'\<oii)i\wi,\ + \vi,\). 



We develop 

\{fM = \{flxU<0{l)(^\vi\ + \V2\ + \pi\ + \p2\ + \vix\^ + 0{l)\viaix\. 

Since 

0'jLO \wi\ = \Vi:, - XiVi + pix - Xipi\ < 0{l)\vi\ ^ \vix\ < 0{l)\vi\ + \pix - Xipi\, 



one has 



\Vlx<^lxVl\ 



\0'vivix\ 



< 0{l)\wixVi - wivi^l + 0{l) \^\wix\ + 0{l)\vix\j \pix - Xipi\. 
Using the previous estimates, we get 

\viaix{fia)x\ < ^(l)!^!^:^! - Vi^Wil + 0{l){\vi\ + \wix\) {\V2\ + \pi\ + \p2\) 

+ 0{l)\wixvi - wivixl + 0{1) (\wix\ + 0{l)\v,x\) \pix - Ami + 0{l)vfx^lxl-^,/,,l<3S} (^)'- 



6. 



\vi{wix - Xiwi) - wl\ = \viwix - vi^wi + vixWi - Xiviwi - wl\ 

< \viWix - VixWil + \wi{vix - XiVi - Wi)\ 

< \viWix - Via;Wi\ + \wiipia; - XiPi)\. 
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7. 



\wixivix - Xivi) + iwiviai)x + (wix - Xiwi)vix\ 

= \2wix{vix - Xivi - wi) - wi^vix + XiWixVi + 2wixWi 
+ wi^viai + wivixcri + wiViai.^ + Wi^^wij; - XiWiVix\ 
= \2wix{pix - Xipi) + 2wixiwi + aivi) - aiwi^vi 

+ XiWixVi + aiWiVia; + [wixVi - WiVix)0'{wi/vi) - AiWiWi:^! 

< 2\wix{pix - Xipi)\ + 2\wix{wi + aivi)\ 

+ |(Ai - ai){wixVi ~ wivix)\ + (wixVi - wivix)0'{wi/vi)\ 

\wi{wi + aiVi -pix + Xipi)\ = \wi{wi + CTiWi) + wi{pix ~ Xipi)\ 



9. 



\wiaix{fla)x\ < \vi(Jlx{'ria)x\, 

and therefore one comes back to case (5). 
This completes the proof of the estimate . 

A. 3 Appendix to Paragraph El 
A.3.1 Proof of the estimate 

It is convenient to introduce a representation formula for pi, i = 1, 2. To this end, two new 
convolution kernels are needed: let 1^^ °{t, s, x) be the solution of the equation 

It' + Aj /a;' Ixx = 0, 

with boundary and initial data 

/^'°(0, s, a;) = 0, s, 0) = ,5t=„ I^' ''{t, s, L) ee 0. 

Without specifying the explicit expression of , we observe that 

/^•* s, x)ds = J^'°{t, x) 



(see equation (|2.6|l for the definition of J^^ Analogously, let 1^^ ^{t, x) be the solution of the 
equation 

if^"- + X*Ix''' -Ixl"- ^0, 

with boundary and initial data: 

/^'0(0, s, x) = 0, s, 0) = 0, I^'"it, s, L)^St^,. 

By construction, it satisfies 

I^'^^{t, s, x)ds = J^'^^it, x) 



(see equation H2.8|l for the definition of J^^ If i < 1 the fimction pi admits the following 
representation formula: 

pi{t, x) ^ / I^'^"{t, s, x)pi{s,0)ds+ I^'^^{t, s, x)pi{s, L)ds+ 

Jo Jo 



A^i (f - s, X, y) ( {X\ - Xx)pxy ~ Xiypi ) (s, y)dyds, 
Jo ^ ' 
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and hence 



+00 



pix{t, x) = / Ix'^^it, s, x)pi{s, 0)ds + Ix'^^it, s, x)pi{s, L)ds 



t i-L 



A 



"'0 



(t - s, X, y) (^{Xl - Xi)piy ~ XiyPi^ (s, y)dyds. 



From the expression of A'^'i , given by formula l|2.4(l . it follows that 



A^(t, •,y)exp(c(. -y)/2) 



and from the previous observations 



exp(ca;/2) j Ix^ ^{t, s, x)ds 


-foo 

Ix' (t, s, x)ds 





tA* 



exp(ca;/2)j;* "(t, a;) < C(l) 



I ji"' "''(t, a;) I 



Hence 

I exp(ca::/2)pi(t, x)\ 



+00 



exp(ca::/2) / Ix^ (t, s, x)pi{s, 0)ds 
Jo 

t 



+00 



exp(ca;/2) / Ix^'"{t, s, x)pi{s, L)ds 



exp(ca;/2) / / Ax^ {t - s, x, y)(^{Xl ~ Xi)piy - XiyPi^{s, y)dyds 



< O{l)\p{x^0)\oo+O{l)\p{x^L)U 



+ 0(1)61 
and therefore 

The estimate 



\ V 



/ JO 



( suppij;(s, y) exp(cy/2) ) / A^ (t - s, x, y) exp (c(x - y)/2)dsdy 



0{l)5l 



I s\xppix{t, x) exp(cx/2)| < 0(1)51 Vi < 1. 

X 

sup|p2x(i, a;)exp (c(L- 2:)/2)| < 0{l)5i Vt < 1 



follows by symmetry. 

If t > 1 the following representation formula holds: 



pix{t,x)= / piit~l, ,y)A^^{l,x,y)dy+ / Ix' \l, s, x)pi{s, 0)ds 
Jo Jt-i 



+00 



tKO, 



-00 



+ Ix^ (1, s, x)pi{s, L)ds 

Jt-i 

Jo Jo 



It follows that 



Ax^ (1 - s, X, y) {^{Xl - Xi)piy - Xiyp^ (t - 1 + s, y)dyd& 
sn^pix{t, x)exp(ca;/2)| < 0{l)6i > 1, 



and hence by symmetry 



I supp2i:(i, x) exp (c(i - x)/2) \ < 0{l)Si Vi > 1. 

X 

This concludes the proof of 1)4.9(1 . 
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A.3.2 Proof of Proposition lOl 

We will perform the computations only for ?J2, W2 and W2x, since those for vi, wi and wix follow 
by symmetry. 

Three new convolution kernels: the solution of equation 

Qt + X*2Qx-Qccx^0 (A.ll) 

with boundary conditions 

Q{0,x)=Sy, g(i, 0) = 0, Qx{t,L)=0 

is 

Q{t, x) = e^'^it, X, y) := (t){t, z, y)(^^G^{z + 2mL - y)+G^{z + 2mL + y)^dz (A.12) 
As in Section 13 we use the notation 

(l)^Ht,x,y)=exp(4ri^-y) ^^'^ ' 



2 ^ 4 
and G{t, x) = exp(— a;^/4t)/2\/7rf. Note that, by construction, 

e;^^'(t, 0, y) = o vt>o, yelo, l[. (a.is) 

Moreover, an argument similar to that used in Section f4 . 1 . II ensures that a maximum principle holds 
for equation (|A.lip . in other words if 

Q(0, x) < 0, Q{t, 0) < 0, Qx{t, L) < 0, 



then Q{t, x) <Oyt, x. 

The solution of (|A.11(I with boundary conditions 



Q(0, x) = 0, Q{t, 0) - 1, QS, L) = 



IS 



B^^{t, x) = l- / e^-it, X, y)dy. (A.14) 





In the following, we will need another convolution kernel, Q^'^{t^ x, y), such that 

QiHt,x,y)^~Q^Jit,x,y). (A.15) 
We arbitrarily impose G)^^ (^t, x, L) = Vt, x and define 



Recalling (|A.13|I . we observe that Q'^'^ {t, x, y) is the derivative with respect to a; of a function z such 
that 

f < X < y 

Zxit, 0, y) = Zx(t, L, ii) = z(0, x, ii) = i . ~ j- 

x\ , , jj x\ , , J) \ 1 1 jj [ 1 y < X < L (A. 16) 

+ -^2^a; ~ Zxx — 0. 

It follows that 8^2 (^t, X, y) satisfies 

e^=(t, 0, y) = e^^(i, L, y) = e^H0,x,y) = 5y 
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and hence actually 

e^Ht,x,y) = A^Ht,x,y), (A.17) 

where A'^^ is the convolution kernel defined by (|2.4|l . In the following, however, for sake of clearness 
we will write Q^^ {t, x, y) when we want to underline that the relation (|A.15|I holds. From the 
identity (jA.lTp and the estimates (|2.11|) it follows 

\\Q^Ht.^.y)\W<0{l) ||e^^(t, X, <^ Vt<l. (A.18) 

Moreover, let z be as in l|A.16|l and let B^'^ be defined by (|A.14|I . then z{t, x, 0) + B^^t, x) = I 
and hence 

e^Ht,x,0) + B^^{t,x)=0. (A.19) 
Such an identity, together with (jA.18|) . implies 

\\B^^t,x)\\Li<0{l) \\B^i{t,x)\\Li<^ t<l. (A.20) 

Since the kernels introduced so far will be used to prove the integrability of V2x with respect to time, 
one has to prove that they are integrable on small time intervals. 

\ei^t, X, y)\dt ^ \A^^t, X, y)\dt < 0(1) Va; e [0, L], Vye]0,L[ (A.21) 
Jo 

Proof. One can check that 

\G^^t, X - y)\dt < 0(1) / \G^'^{t, x-y)\dt<Oil) Vx, y G M. (A.22) 



^0 



Since 

, a; 



(i, x,y)=[ (b{t, X, y) ^ G{t, x-y + 2mL) \ - f x, y) ^ G{t, x + y + 2mL) 

m>0 ' ^ m>0 

<^(t, X, y) G{t, x-y- 2nL) ) - U(t, x, y) G{t, x + y- 2nL) 



n>0 ^ ^ n>0 

one gets 



{t, X, y)\ < |G7 {t,x~y + 2mL)\ + ^ |G^ {t,x + y + 2mL)\ 

m>0 m>0 

+ Y\G^-Ht, 2nL + y-x)\+X;Y\G^Ht, 2nL + y - x)\ 

n>0 n>0 

+ Y\G^-Ht, 2nL-2/-x)|+A2^|G^^*(t, 2nL - y ~ x)\. 



n>0 n>0 

Since 

\G^Ht, z + 2mL)\ < e-'"-^|G'^'(i, z)\ \G^Ht, z + 2mL)\ < e-"'^\G^^t, z)\ 

if m > 0, i < 1 and z is large enough, from the previous estimates and from (|A.22|I one deduces 
(TOTll. □ 



From equation ljA.19p and the previous estimate it follows 

x)\dt<0{l) Vxe[0, L]. (A.23) 



Jo 
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A representation formula for V2 : it is convenient to introduce the auxiliary function 



which satisfies the equation 
where 



V2{t, X) = / V2{t, 

Jo 



1^24 + -1^2.. =5i(i, X), 



SAt, x) = / hit, OdC 
lo 



The boundary and initial conditions of V2{t, x) are 



¥2(0, X)= «2(0, V2it, 0) = iv2., - X2V2){S, 0)ds, V2.(<, i)-0 

Jo Jo 
The convolution kernels ljA.12|) and ljA.14|) provide the representation formula 

V2{t, x)= I (i, X, y)V2{Q, y)dy + [ B{t~ s, x){v2^ - A2i^2)(s, 0)ds 
Jo Jo 

{t - s, X, y) ( (a; - A2)w2 ) (s, y)dyas 

lo Jo 



{t - s, X, y)(^{X*2 - A2)w2) (s, y)dyds 
_ s, X, y)Si (s, y)dyds. 



(A.24) 



Jo 



Since 

from l|A.24p one gets 



e^^ {t, X, 0) + B^^ (t, x) = 0, Si {t, 0) = 0, 



V2At, ^V2{t, x) = / e^^( 

Jo 

e^'^{t- s, X, y)si{s, y)dyds + / e^-^(t - s, x, y)\[A2 - ^2)V2 } (s, y}ayas 

Jo Jo ^ ^ y 



'■{t,x, y)v2 (0, y)dy+ j B^' {t~s, x) (v2:, - \*2V2){s, 0)ds 

Q^'^{t-s, X, ?/)((A^ - A2)f2)^(s, y)dyds 



and 



V2x(t, x) ^ Qx'^{t, x, y)v2{Q, y)dy + I B^l{t - s, x)[v2x - Kv2]{s, Q)ds 
Jo Jo 

ne^" [t - s, X, y)si(s, y)d2/ds 



e;^^(i - s, X, y)(^{X2 - \2)v2y " A2y'y2) (s, y)dyds 



From the estimate (lA.lSp , (jA.21|) and ljA.23|) on the convolution kernels it follows 

1 .1 
\v2xit, x)\dt < |lf2(0)|Ui sup / |e^2(<, X, y)\dt 

x,V Jo 

+ ( I \{v2x - \2V2){s, 0) + (A; - A2)i;2(s, Q)ds\ds 



< 0{l)6i 
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for all X €z [0, L]. If t > 1 we can use for V2x the expression 

V2xit,x)^ f e^^l,x,y)v2{t-l,y)dy+ f 3^1(1 ^ s, x){v2, - X*2V2){t ~ 1 + s, 0)ds 
Jo Jo 

"1 rL 

e^' (1 - s, X, y)h{t -l + s, y)dyds (A.25) 



'0 Jo 

Computations analogous to the previous ones lead to 

\v2As, x)\ds < 0{l)6i. 

Ji 
Hence 

i-T 



\v2x{s, x)\ds <0{l)5i VT>0, xe[0,L]. 
The integrability of W2 with respect to time: it holds 

\w2{s, y)\ds < 0{l)di yt>0, yye[0,L]. (A.26) 



Proof. We preliminary observe that 

W2{0, X) = {£2, Ut{0, X)), ■W2{t, 0) ^ {I2, u'f,o{t)), W2{t, L) ^ {£2, Ui,j^{t)), 

where £2 satisfies {£2, 7'2) = 1 and {£2, ri) = 0. Hence 

\\w2{t = 0)||l1(o,L) < 0{l)Si, \\W2{X = 0)||l1(0,+oo) < Si, \\W2{X = i)||Ll(o,+oo) < ^1- 

Let 2c be the separation speed defined by (|1.5|) , let if be a compact neighborhood of the value u* 
defined by H1.8|l and let C > satisfy 

0<c<A2(u)<C yueK. 

If y G]0, L[, the estimate (|A.26ll can be obtained applying Lemma [4. II to the functional 

{a(l - 6^"-"=^) X < y 
(A.27) 

where a and b satisfy 

a(l - e-'^y) = b^e-^y - e-"^) 

(A.28) 

aCe'^y + bce-''y = 1. 
By straightforward computations, from ljA.28(l one gets that the functional Py satisfies 

Py{0)^Py{L)^0, Q<Py{x)<Py{y)<0{\), P'y{Q)<0{\), -P'y{L) < 0{\), VL » 1 

P'^{x)+\2P'y{x)<-6x^y. 

Since W2 satisfies 

W2t + {Mw2)x - W2xx = hit, x), 

Lemma [4.11 ensures that 

\w2{s, y)\ds < 0{1) \w2{0, x)\dx + 0{1) / {his, x)\dxds 



Jo Jo 

t pt 



+ 0{l)j \w2{s,Q)\ds + 0{l) I \w2{s,L)\ds 
<0{l)6i V2/e]0, L[. 



□ 



49 



Integrability of W2x with respect to time: it holds 



\w2xis, x)\ds <0{l)Si Vt>0. 



(A.29) 



Proof. From the representation 



t rL 



W2x{t, x) ^ I Ap {t, X, 2/)w2(0, y)dy +11 A^" {t - s, x, y)s2(s, y)dyds 
Jq Jo Jo 

Ax{t - s, X, y)\{\*2 - ^2)'W2y - X2yW2){s, y)dsdy + W2(0, L)Jx^ {t, x) 
'0 Jo ^ ' 

+ W2(0, 0)J^"(t, a;) + / Jx'°it- s, x)w'2{s,0)ds+ J^' ^ (t - s, x)w'^{s, L)ds 



it follows 

If t > 1 one can write 



\w2x\{t, x)dx < 0{l)Si 



W2xit, x) ^ I Ax^{l, x,y)w2it- I, y)dy+ I I A^^ {1 - s, x, y)s2it - I + s, y)dyds 
Jo Jo Jo 

+ J J Ax^{l ~ s, X, y)(^{\2- \2)w2y - X2yW2^{t - 1 + s, y)dsdy 

+ W2{t - 1, L)Jx^ ^(1, X) + W2{t - 1, Q)Jx'\l, X) 

+ I Jx" °(1 - s, x)w'2{t - 1 + s, 0)ds + / Jx^ ^(1 - s, x)w2{t - 1 + s, L)ds 



1 



and obtains 



\w2x\{t, x)dt<0{l)5i 



This concludes the proof of ljA.29|l . 

A. 3. 3 Proof of the estimate 

We need three preliminary results: 
• For any t <1, the following holds: 



□ 



\Qf{t,x,y)\<a[t,x^y) + b{t,x) \\a{t)\\m-L.L). l) < (A.30) 

Proof of (|A.30|I In the following, a(t, x — y) and x) will denote functions that satisfy 

0{1) 



0{l) 



By the identities TT^i and HA.17|) . 

exHt,x,y)^^xHt,x,y)=[<l>^Ht,x,y) ^ G(i, x + 2™^ - y) - G(i, x + 2mi + y) 
One has 

m — +00 



m — +00 



m — — 00 



m — — cxD 



y) ^ G{t,x + 2mL-y)j < ^ | (i, x - y + 2mL) + A2 ^ G^^ (t, 2nL - x + y) 

m>0 n>0 

^ Gx^ {t, 2nL ~x + y) <a{x~y), 



n>0 
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where we have set n := — m. To complete the proof of (|A.30p . it is convenient to observe that 

G^^ (t, X + y) < G^'^ (t, x) yx> {Xlt + V2t}, Vy > 

and that 

\G^Ht,^ + y)\<GtHt, x) + GAt, ^^)X|o< :c<V2i + x;t}-^^'^^ 

IG^^t, 2L-x~v)\<G^-[t,L-x) + G.,{t, ^2^) X| ^ _ ^ _ < ^ < < /3(a:), Vx, y e [0, 
where xe denotes the characteristic function of the set E. Hence 

m — + 00 N 

^'^{t,x,y) G{t,x + 2mL + y)\ < ^ G^^ (i, a; + y + 2mL) + G^' (t, x + y) 



m — — oo 



/ X 



m>0 



+ ^ G^2 {t^ 2ni - X - y) + ^ G^ [t, 2nL - x - y) 

n>0 n>0 

< XI ^ + 2^"^) + + A2 XI '^^^ i - a;) + X (t, (2n - 1)L - a;) + Gx^ (i, 2L - 

m>0 ri>0 n>l 

which concludes the proof of l)A.30|l . □ 
• If t < 1 then 

x)\dx < (A.31) 
vt 

Proof. Let t < 1. From the equality 

Uxx = Vi(l)fiUx +Vixfiy +(7lxfla^ + Vlxh + Pl(BfiUx + pi^fip^ + Plxh + V2xr2 +P2x?'2, (A. 32) 

and from the bounds ||pi2;(i)||Li < C(l)^i and < 0{l)5i/ ^/t, it follows that 

\\vix{t)\\ = 11(4, U.,x(t)) -pl.WIUi < 

where £i = (1, 0). Hence 

\\w,mL^ < 0(l)||«l(t)|lL^ + \\vix{t)\\L^ + 0{l)\\p,(t)\W + \\Plx{t)\W < 



From the estimates 



\\w'l{x = 0)||l1(0,+oo) 11(4, MbO>IUi(0,+oo) < h 
\\w'i{x = i)||Li(0,+oo) = 11(4, UhL)||Li(0,+oo) < -^l 

= 0)||li(o,l) = 11(4, Wo-^(«oK)||lho,l) <^^(l)'5i, 



and from the representation formula 



x) = / A^i(t, x, y)u;i(0, y)dy + J^''°{t, x)wi(0, 0) + J,^*^(t, a;)u;i(0, L) 





+ / Jx'°{t- s, x)w[{s,0)ds+ [ Jx'^it- s, x)w[{s, L)ds (A.33) 
Jo Jo 

A:r'(*- s, a;, y)((Ai - Ai)wiy - XiyWi){s, y)dsdy, 
Jo ^ ^ 
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it follows that 



Hence 



ki,(i)llLi < 



and therefore from (|A.32|I one gets (|A.31 



Vl 



< 



Oil)Si 

Vt 



• If <> 1 then 



□ 



(A.34) 



Proof. One can repeat the same computations performed to prove IjA.Sip . using, instead of (jA.33|l . 
the following representation formula (which holds if t > 1): 



Wlx{t, x) 



At' (1, X, y)wi{t - 1, y)dy + J^^ "(1, x)wi{t - 1, 0) + J^^ ^(1, x)wi{t - 1, L) 



+ / Jx'^ "(1 - s, x)w[{t - 1 + s, 0)ds + / Jx' ^(1 - s, x)w[{t - 1 + s, L)ds 
Jo Jo 

+ Ax'{l - s, X, y)(^{Xl - Xi)wiy - XiyWi^{t - 1 + s, y)dsdy. 



□ 



Let 



X(r) := sup 

re(-T,T) 



min{T, T+r} ^min{L,L+^} 



<{0,t} 



<{o,a 



|w22;(^ — T, x — ^)\dtdx. 



dt < 0{l)Sl 



It holds: ^ ^ 

/ / \vi{t,x)\\v2xit,x)\dxdt<IiT). 
Jo Jo 

Moreover, thanks to the estimates (|A.31|I and (|A.34p . 

r.max{2, 2+t} (.min{L, 
<{0, r} Jmax{0,?} 

Hence we are left to estimate the term 

pmin{T,T+T} pmin{L,L+i} 
Jmax{2,2+T} Jmax{0, 

in the case T > 2: to do this, we will exploit the representation formula (|A.25ll and the estimate 



\vi(t, x)\\v2x{t-T,X-0\<O{l)\\v,\\L^.di ^'|l + i=| 



\vi{t, x)\ \v2x{t -T, X - ^)\dxdt 
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One has 

min{T. T+r} ^min{L,L+?} rL 

max{2,2+T} Jmax{0, J} JQ 

<.min{T, T+t} <.min{L, <.L 

< / / vi{t,x)l a{l, X - ^ - y)v2{t - 1 - T, y) 

v'max{2,2+T} Jmax{0,{} JO 

inin{T, T+t} <.min{L, L+C} /-L 

/ vi{t, x) / 6(1, x-e)i'2(t-l-r, ?/) 

max{2,2+r} Jmax{0, ^} JO 
L <.min{T, T+t} <.min{L, L+J+z} 

< / a(l, / / vi{t, x)v2it - 1 - X - z - £,)d^ 

L Jmax{2,2+r} Jmax{0, 2+^} 

min{L,L+5} / <.min{T, T+r} / /"^ \ \ 

b{l,x-0[ / t-ili, 2:) / f2(t-l-r, da;<0(l),S2^ 

max{0,4} YJmax{2,2+T} \J0 J J 

and 

^min{T, t+t} <.min{L, L+5} pi pL ^ 

I vi(t, x) / 6a;^ (1 - s, x - ^, ?;)( (A2 - A2)w2y) (t - T - 1 + s, ?/)dydsda;d< 

max{2,2+T} Jmax{0, C} JO JO ^ ' 

^inin{T, T+t} ^min{L, L+5} ^1 pL 

< 61 / tii(t, x) / / a{l — s, X — — y)v2y{t — T — 1 + s, y)dydsdxdt 

Jmax{2,2+T} Jmax{0, 5} JO 

<.min{T, T+t} pniin{L,L+^} pi pL 

+ Si / vi{t, x) / / b{l — s, X — S,)v2y{t — T — 1 + s, y)dydsdxdt 

Jmax{2,2+T} Jmax{0, 5} Jq 
pi pL I <.min{L, L+z+5} ^inin{T, T+t} \ 

< 61 / a(l — s, z) / / vi{t, x)v2x{t ^ T — 1 + s, X ~ £^ — z)dxdt \dzds 

JO J-L Wmax{0,C+2} Jmax{2, 2+t} / 



<5i 



nmin{L,L+5} / <.inin{T, T+t} / /"^ \ \ 

5(1 — s,a; — ^) / vi[t^ x){ I V2y{t — t ~ 1 + y)dy\dt\dxds 

jax{0, 5} \ Jmax{2, 2+t} \J0 ) j 

<0(l)<5iX(T)+0(l)J^ 

With analogous computations one can estimate the other terms that comes from the representation 
formula ljA.25p and hence prove that I(T) < 0(\)b\. 

A.3.4 Proof of the estimate (|IT7| 

Since in the following we will often refer to equations H3.1()|l and 13.1111 , we recall them: 
where the cut-off d is given by 

r s if Ml < -5 

Q{s)^l if|s|>3<5 h«l<-. 
\ smooth connection if (5 < s < 3i5 

It follows that \w\ + a\V\\ ^ only when the function 9 is not the identity, i.e. when \wi + A*wi| > 
5\vi\. Since 

wi = vi^ - Xivi + pix - XiPi, 
the condition \wi + aiVi \ implies 

\vix + (Ai - Ai)vi + pix - Aipil > S\vi\. 
There are therefore two possible cases: 
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1. 

\vix + {K - Ai)t.i| > -5\vil 
and therefore, since |A* — Ai| < 0{l)5i and 5i « S, 

\vix\ > ^\vi\. 

2. 

S 6 

\vix\ < -^\vi\ =4> bla; - Aipil > 

If case 1 holds, then 

\wi + aivi\ = \vix + ((Ti - Ai)ui +pix - Aipil 

< \vix\ + Sl\vi\ + \pij; - Aipil < 0{l)\vix\ + Ipi-j: - Aipil 

and therefore 

\vi\ + \wi\ + \vi^,\ + \wi \){\w,+a^v,\) <0(1)(| Vlx\ + \pi\ + \vix\ + \wix\] {0{l)\vix\ + \pix - AlPil 
< 0{l)(\vi.,\ + \pi\ + \pix\ + ki.l) \pix - AlPil + + (bil + + 

Since 

bil, \vxx\ < 0(l)^iexp(-ca;/2), 
it follows that, if case 1 holds, then one is left to prove 



X {|(.,M)+Atl>5} + l^^i-l") ^)d^dt < 0{l)6l (A.35) 



/O ^0 

On the other hand, if case 2 holds then 

4 - 

\vix + (ci - Ai)wi + pix - Xipi\ < -5\vi\ + \pix - Xipi\ < 0(1)1^1^; - Aipil, 
and therefore 

J J (\vi\ + \wi\ + \vix\ + \wia;\^(^\wi+aivi\^{s, x)dsdx 

<0(1)^ ^ {\vi\ + \wi\ + \vix\ + \wix\^\pi,-XiPi\is,x)dsdx<0{l)Sl 
thanks to the exponential decay of \pi \ and \pix \ ■ 

To prove l|A.35|) it is convenient to introduce a new cutoff function: 

r if \s\ < 3/5(5 
ip{s)^l 1 if |s|> 4/5(5 

I, smooth connection if 3/5 S < \s\ < 4/5(5. 

Moreover, in the following we will only prove that 

^'^^''x{|(.,/.,)+Atl>^}l^i-l'(^'^)^^'^^^'^(l)'^?' (A-36) 

because the estimate 

X{\(n^,,/v,)+xi\>s}\wixf{t,x)dxdt < 0{l)Sl 



"'0 
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can be obtained with similar techniques. 

As we have already observed, it is sufficient to show 



T i-L 



^0 



Multiplying the equation 



+ At) (t, x)dxdt < 0{l)Sl 



Vit + (AiWi)a; - Vixx = 



by ipvi, we get 




d ,vl 



Jo 

T 



{-^^) ^(V't + ^l'<Px - i>xx) + ^Vix\^ + -Y^lxlp - vli)xx 1 dxdt 



dt^ 2 







i)Vi{\iVi - Vix) 



x—L 



dt 



x=0 



-j-{^x~\l^) 



X—L 



(A.37) 



dt. 



x=0 



Indeed, 



and 



d v'^ v'^ 



L nT 



L rT 



(Aiui - vix) ^i^vidxdt ^ / (vix - Xivi){-)pvi)xdxdt + 



-1 x—L 



^Vi{XiVi - Vix) 



dt 



L i-T 



'0 JO 



^x\^] +iJvl^-\ii}xvl-Xii}\^ \ dxdt + 



tpvi{XiVi 



dt 



Jo 



V'wL + ( ^ ) i^ixip - Xiipx + tpxx - 2-4)xx)dxdt + 



x=0 
-1 x—L 



^l}Vi{XiVi - Vix) 



dt 



x=0 



-jii'x - Xi^l^) 



dt. 



x=0 



One can develop the term V't + Xitpx — i'xx and, since 



[ vl 



fpx = Ip' 



I WlxVl - WiVix 



Wi 



I ( W \xxV\ - VlxxWl _ ^ VlxjwixVl - WiVix) 

n ^ n 



(A.38) 



'<Pxx = ^y—j +^ 
one obtains 

vl{%pt + Xilpx - Ipxx) = 1p'vi{wit + {XiWi)x - Wixx) - Ip'wiivit + {XiVi)x - Vixx) 



III 2r^i 



2i}j'vixVi 



Thus, inserting the last formula into ljA.37|l . we obtain 



L I-T 



^0 



L r 1 *-=T 

2, 



Vidx 

L I-T 



t=o Jo 



'4'Vi{vix - XiVi) 



X—L 



dt 



-I X — L 



^(V'.-Ai^) 



dt 



III 2 1 ^1 \ // / '"1 \ 2 / '^1 \ / 

Vj — +ipvixVi[ — + v-f^tpxx - -;:rXixV 



Wi 
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The boundary terms are bounded by 0{l)5f since HwiUl^o < 0{l)Si and thanks to the estimates of 
Proposition Ol Since by (|4.16|) 



we are left to estimate the fohowing terms: 



T pL 



,1 I 



dsdx < 



T pL 



,1 I wi 

W Vlx Wlx Vlx 

Vl 



-0' ( \vi\ + \pix - Aipi I ) ( Wix - — 



dsdx 
dsdx 



< Oil) I I 

JQ JO 

<0{1) [ [ \viwix - vixwi\dsdx + 0{1) [ [ \pix ~ XiPi\(\wix\ + 0{l)\vix 

jQ Jo Jo JO ^ 

Indeed, if ip' ^ then |A* ~ wi/vi\ < 6 and hence 



T j-L 

Jo 



'^1^' [wixxVi - wivixx^dsdx = 

<0(1) 

L I tHHX i >-A;w- 
(A,-A;)(-|* 





■ 




r 


lo J 









f 


Wi 


Jo 





(^wixVi — wivix^ dsdx 



WixVi - WiVix 



X—L 



<0{l)Sl 



x=0 



'■L rT „2 



L i-T 



x=0 



2 ^ 



+ Oil)6f< 0(1)61 



x=0 
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